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We prove that the system of subordination equations, defining the free additive con¬ 
volution of two probability measures, is stable away from the edges of the support 
and blow-up singularities by showing that the recent smoothness condition of Kargin 
is always satisfied. As an application, we consider the local spectral statistics of the 
random matrix ensemble A + UBU*, where U is a Haar distributed random unitary 
or orthogonal matrix, and A and B are deterministic matrices. In the bulk regime, 
we prove that the empirical spectral distribution of A UBU* concentrates around 
the free additive convolution of the spectral distributions of A and B on scales down 
to Af-2/3. 
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1. Introduction 

One of the basic concepts of free probability theory is the free additive convolution of two 
probability laws in a non-commutative probability space; it describes the law of the sum of 
two free random variables. In the case of a bounded self-adjoint random variable, its law 
can be identified with a probability measure of compact support on the real line. Hence 
the free additive convolution of two probability measures is a well-defined concept and it is 
characteristically different from the classical convolution. 

In this paper, we prove a local stability result of the free additive convolution. A direct 
consequence is the continuity of the free additive convolution in a much stronger topology 
than established earlier by Bercovici and Voiculescu [10]. A second application of our sta¬ 
bility result is to establish a local law on a very small scale for the eigenvalue density of a 
random matrix ensemble A -\- UBU* where C/ is a Haar distributed unitary or orthogonal 
matrix and A, B are deterministic A by A hermitian matrices. 

The free additive convolution was originally introduced by Voiculescu [36] for the sum 
of free bounded noncommutative random variables in an algebraic setup (see Maassen [32] 
and by Bercovici and Voiculescu [10] for extensions to the unbounded case). The Stieltjes 
transform of the free additive convolution is related to the Cauchy-Stieltjes transforms of the 
original measures by an elegant analytic change of variables. This subordination phenomenon 
was first observed by Voiculescu [38] in a generic situation and extended to full generality by 
Biane [14]. In fact, the subordination equations, see (2.5)-(2.6) below, may directly be used 
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to define the free additive convolution. This analytic definition was given independently by 
Belinschi and Bercovici [4] and by Chistyakov and Gdtze [18]; for further details we refer to, 
e.g., [39, 27, 2], 

Kargin [30] pointed out that the analytic approach to the subordination equations, in 
contrast to the algebraic one, allows one to effectively study how free additive convolution 
is affected by small perturbations; this is especially useful to treat various error terms in the 
random matrix problem [31]. The basic tool is a local stability analysis of the subordination 
equations. In [30], Kargin assumed a lower bound on the imaginary part of the subordination 
functions and a certain non-degeneracy condition on the Jacobian that holds for generic 
values of the spectral parameter. While these so-called smoothness conditions hold in many 
examples, a general characterization was lacking. Our first result. Theorem 2.5, shows 
that the smoothness conditions hold wherever the absolutely continuous part of the free 
convolution measure is finite and nonzero. In particular, local stability holds unconditionally 
(Corollary 2.6) and, following Kargin’s argument [30], we immediately obtain the continuity 
of the free additive convolution in a stronger sense; see Theorem 2.7. 

The random matrix application of this stability result, however, goes well beyond Kargin’s 
analysis [31] since our proof is valid on a much smaller scale. To explain the new elements, 
we recall how free probability connects to random matrices. 

The following fundamental observation was made by Voiculescu [37] (later extended by 
Dykema [20] and Speicher [35]): \i A = and B = are two sequences of Hermitian 
matrices that are asymptotically free with eigenvalue distributions converging to probability 
measures and fip, then the eigenvalue density oi A + B is asymptotically given by the 
free additive convolution fj,a H iJ-p- One of the most natural ways to ensure asymptotic 
freeness is to consider conjugation by independent unitary matrices. Indeed, if A and B 
are deterministic (may even be chosen diagonal) with limit laws fXa and then A and 
UBU* are asymptotically free ii U = is a Haar distributed matrix; see [37] and many 
subsequent works, e.g., [34, 40, 15, 33, 19]. In particular, the limiting spectral density of 
the eigenvalues oi H = A + UBU* is given by /Tq, H M/3- 

The conventional setup of free probability operates with moment calculations. An alter¬ 
native approach [33] proves the convergence of the resolvent at any fixed spectral parame¬ 
ter z G C"*". Both approaches give rise to weak convergence of measures, in particular they 
identify the limiting spectral density on macroscopic scale. 

Armed with these macroscopic results, it is natural to ask for a local law, i.e., for the 
smallest possible (iV-dependent) scale so that the local eigenvalue density on that scale still 
converges as N tends to infinity. Local laws have been somewhat outside of the focus of free 
probability before Kargin’s recent works. After having improved a concentration result for 
the Haar measure by Chatterjee [17] by using the Gromov-Mihnan concentration inequality, 
Kargin obtained a local law for the ensemble H = A + UBU* on scale rj ^ (log [29], 
i.e., slightly below the macroscopic scale. Recently in [31], he improved this result down 
to scale rj ^ under the above mentioned smoothness condition. In Theorem 2.8 we 

prove the local law down to scale rj = Imz ^ any additional assumption. 

To achieve this short scale, we effectively use the positivity of the imaginary parts of the 
subordination functions by localizing the Gromov-Milman concentration inequality within 
the spectrum. Since the subordination functions are obtained as the solution of a system 
of self-consistent equations whose derivation itself requires bounds on the subordination 
functions, the reasoning seems circular. We break this circularity by a continuity argument 
(similarly as in [23]) in which we reduce the imaginary part of the spectral parameter in 
very small steps, use the previous step as an a priori bound and show that the bound does 
not deteriorate by using the local stability result. Theorem 2.6. 
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Finally, we remark that the local stability result is also a key ingredient in [3], where we 
were able to prove a local law down to the smallest possible scale rj ^ N~^, but with a 
weaker error bound than in Theorem 2.8; see Remark 2.4 for details. 

1.1. Notation. We use the symbols O(-) and o( •) for the standard big-0 and little-o 
notation. We use c and C to denote positive numerical constants. Their values may change 
from line to line. For a,b > 0, we write a < 6, a > & if there is C > 1 such that a < Cb, 
a > C~^b respectively. We write a ^ 6, if a < & and a>b both hold. We denote by ||r!|| the 
Euclidean norm of £ C^. For a.n NxN matrix A £ Mn{C), we denote by ||^|| its operator 
norm and by ||A ||2 := \/(A, A) its Hilbert-Schmidt norm, where {A,B) := Trace(kli?*), for 
A,B G Finally, we denote by tr^ the normalized trace of A, i.e., trA = ^Trace^. 

Acknowledgment. We thank an anonymous referee for many useful comments and re¬ 
marks, and bringing references [7, 12] to our attention. 


2. Main results 


2.1. Fbee additive convolution. In this subsection, we recall the definition of the free 
additive convolution. Given a probability measure* /i on K, its Stieltjes transform, on 
the complex upper half-plane C"*" := {z £ C : Im^ > 0} is defined by 


.W := / 


d^(a;) 

X — z 


£ C+. 


We denote by the negative reciprocal Stieltjes transform of p, i.e., 
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Observe that 


lim 






z£C+. 


77/^00 IT] 


= 1 . 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


as follows easily from (2.1). Note, moreover, that is an analytic function on C"*" with 
non-negative imaginary part. Conversely, if F : C"*" —>■ C'*' is an analytic function such that 
lim,j/>too F(iri)/ir] = 1, then F is the negative reciprocal Stieltjes transform of a probability 
measure p, i.e., F{z) = F^{z), for all z £ C+; see, e.g., [1]. 

The free additive convolution is the binary operation on probability measures on R char¬ 
acterized by the following result. 


Proposition 2.1 (Theorem 4.1 in [4], Theorem 2.1 in [18]). Given two probability measures 
Pi and p 2 on R, there exist unique analytic functions, wi, W2 : C'*' —>■ C'*', such that, 

{i) for all z £ C’*', Ima;i( 2 :), lmui 2 {z) > Imz, and 


lim 

77/^00 


^i(W) 

ip 


lim 

Tj/^OO 


0 ^ 2 ( 11 ?) 

ip 


= 1 ; 


(a) for all z £ C’*', 


Ffj,^{uJ 2 {z)) - UJi{z) - W 2 {z) + z = 0, 
F), 2 (u;i(z)) - wi(z) - uj 2 iz) + z = 0. 


(2.4) 


(2.5) 


*A11 probability measures considered will be assumed to be Borel. 
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It follows from (2.4) that the analytic function F : C“^ —^ C“^ defined by 

P{z) ■= = Ff,^{uJi{z)), ( 2 . 6 ) 

satisfies (2.3). Thus F is the negative reciprocal Stieltjes transform of a probability mea¬ 
sure called the free additive convolution of jii and ^ 2 , usually denoted by /i = ffl /i 2 . 
Note that (2.6) shows that the roles of /ii and ^2 are symmetric and thus /ii ffl/r 2 = /r 2 ffl/ii. 
The functions uti and UJ 2 of Proposition 2.1 are called subordination functions and F is said 
to be subordinated to , respectively to F)j 2 . 

We mention that Voiculescu [36] originally introduced the free additive convolution in 
a different, algebraic manner. The equivalent analytic definition based on the existence of 
subordination functions (taken up in Proposition 2.1 above) was introduced in [4, 18]. 

We next recall some basic examples. Choosing fj,i arbitrary and fi 2 as a single point 
mass at 6 G M, it is easy to check that fii ffl ^2 simply is fXi shifted by b. We exclude this 
uninteresting case by henceforth assuming that and /i 2 are both supported at more than 
one point. Choosing /ri = ^2 = A* as the Bernoulli distribution 

At= (1 - ^),5o-f f G (0,1), 

the free additive convolution is explicitly given by (see e.g., (5.5) of [33]) 

(/i ffl ^,){x) = y(^+ - ^)+(^ - { zl± + (1 _ 20Mx) + {2f - 1)+<52(x) , (2.7) 

ivxyl — x) 

a; G K, where := 1 ± — ^) and where (•) + denotes the positive part. Observe that 

/r ffl ^ has a nonzero absolutely continuous part and, depending on the choice of a point 
mass. Another important choice for fj ,2 is Wigner’s semicircle law Hsc- For arbitrary fii, 
fii ffl fisc is then purely absolutely continuous with a bounded density^ that is real analytic 
wherever positive [13]. 

Returning to the generic setting, the atoms of ffl ^2 are identified as follows. A point 
c G M is an atom of fii ffl /r 2 , if and only if there exist a,5 G M such that c = a + b and 
Mi({a}) + M 2 ({^}) > 1; see [Theorem 7.4, [11]]. For another interesting properties of the 
atoms of fii ffl fi 2 we refer the reader to [12]. The boundary behavior of the functions 
^2 has been studied by Belinschi [5, 6, 7] who proved the next two results. 
For simplicity, we restrict the discussion to compactly supported probability measures. 

Proposition 2.2 (Theorem 2.3 in [5], Theorem 3.3 in [6]). Let fii and ^2 be compactly 
supported probability measures on R, none of them being a single point mass. Then the 
functions uii, 1 x 2 ■ C"*" —?► C’*' extend continuously to R. 

Belinschi further showed in Theorem 4.1 in [6] that the singular continuous part of /ii ffl/r 2 
is always zero and that the absolutely continuous part, (/ii ffl of pi ffl fi 2 is always 

nonzero. We denote the density function of {fii ffl fJ. 2 )^^ by 

We are now ready to introduce our notion of regular bulk, of pi ffl pt 2 - Informally, 

we let be the open set on which /ii ffl fi 2 admits a continuous density that is strictly 

positive and bounded from above. For a formal definition we first introduce the set 

■= int i supp {pi ffl ( 12 )^^ \ {x G R : lim + i^) = 0} i . (2.8) 

f v\o j 

Note that does not contain any atoms of pi ffl p, 2 - By Privalov’s theorem the set 

{x G R : lim^x^o F)jjH^ 2 (x -I- ip) = 0} has Lebesgue measure zero. In fact, an even stronger 
statement applies for the case at hand. Belinschi [7] showed that if x G R is such that 


^All densities are with respect to Lebesgue measure on R. 
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lim^x^o Ffj.iSii 2 [x + ir]) = 0 , then it must be of the form x = a + b with //i({a}) +/i 2 ({ 6 }) > 1 , 
a, 6 S R. There could only be finitely many such x, thus must contain an open non¬ 

empty interval. 

Proposition 2.3 (Theorem 3.3 in [ 6 ]). Let fii and fi 2 be as above and fix any x £ 

Then loi, ui 2 : C'*' —>■ C'*' extend analytically around x. In particular, the density 

function fij,,,mij ,2 *'S xeal analytic in wherever positive. 

The regular bulk is obtained from by removing the zeros of inside 

Definition 2.4. The regular bulk of the measure g-i ffl /i 2 is defined as the set 

- \ ^ ■ /miE0M2(®) ■ (2-^) 

Note that is an open non-empty set on which fii ffl ^2 admits the density 

The density is strictly positive and thus (by Proposition 2.3) real analytic on 


2.2. Stability Result. To present our results it is convenient to recast (2.5) in a compact 
form; For generic probability measures/ri, p 2 as above, let the function : (C+)^ — 

be given by 


(W1,W2,^) 


T)ii(a;2) - uji - 0J2 + z \ 
Ffj,^{u}i) - Wi - 0 J 2 + z ) 


( 2 . 10 ) 


Considering as fixed, the equation 


^/il ,/i2 (^1 1 ^2 ; O 5 


( 2 . 11 ) 


is equivalent to (2.5) and, by Proposition 2.1, there are unique analytic functions lv\,u :2 '■ 
C+ —>■ C+, z !->■ uji(z),uj 2 {z) satisfying (2.4) that solve (2.11) in terms of z. We use the 
following conventions: We denote by wi and W 2 generic variables on and we denote, with 
a slight abuse of notation, by uJi(z) and uj 2 {z) the subordination functions solving ( 2 . 11 ) in 
terms of z. When no confusion can arise, we simply write $ for 
We call the system (2.11) linearly S-stable at (a;i,a; 2 ) if 


^iJ.iW2 (‘^ 1 ,^ 2 ) 


-1 

F;ja;i)-1 -1 



( 2 . 12 ) 


for some constant S. Especially, the partial Jacobian matrix, D<l’(a;i, W 2 ), of (2.10) given by 


D$(wi, W2) 


5$ , 

—(wi,u;2, 

UUJi 


z) 


, -—(a;i,W2,2;) 

OUJ 2 


-1 E;^(u;2)-1\ 

F',M)-^ -1 h 


admits a bounded inverse at (wi,a; 2 ). Note that D$(wi,a; 2 ) is independent of z. 

Our first main result shows that the system (2.11) is linearly stable and that the imaginary 
parts of the subordination functions are bounded below in the regular bulk. We require some 
more notation; For a,b > 0, b > a, and an interval I C R, we introduce the domain 


Si{a, b) := {z = E + irj G C“'' : EGT,a<r]<b}. (2-13) 


Theorem 2.5. Let /ii and ^2 be eompaetly supported probability measures on R, and assume 
that neither is supported at a single point and that at least one of them is supported at more 
than two points. Let I C be a compact non-empty interval and fix some 0 < pM < 00 . 

Then there are two constants k > 0 and S < 00 , both depending on the measures fii 
and fi 2 , on the interval X as well as on the constant pM, such that following statements hold. 
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(i) The imaginary parts Imwi and Ima ;2 of the subordination functions associated 
with pii and p 2 satisfy 

min lui(jJi{z) > 2k, min lvauJ 2 {z) > 2k. (2-14) 

2G5x(0,r;M) zeSx(0,r/M) 

(ii) The system = 0 is linearly S-stable at (uji{z),uj 2 {z)) uniformly in 

Sx{0,r]M), i-e., 

max rfj,^^fj,^{uji{z),uj 2 {z)) < S. (2-15) 

2e5i(0,r;M) 

Remark 2.1. The assumption that neither of /ii, /i 2 is a point mass guarantees that the 
free additive convolution is not a simple translate. The case when both, fii and fi 2 are 
combinations of two point masses is special and its discussion is postponed to Section 7. 

Theorem 2.5 has the following local stability result as corollary. 

Corollary 2.6. Let pLi, p ,2 and 5i(0,?7m) be as in Theorem 2.5. Fix zq € C"*". Assume that 
the functions wi, 0 J 2 , ri, ¥2 : C'*' —>■ C satisfy Imd;i( 2 ;o) > 0, l'muj 2 {zo) > 0 and 

(^ 0 ),^ 2 ( 2 : 0 ), Zq) =r{zo), (2.16) 


with r{z) := {ri{z),r 2 {z))~^. Let uJi, u )2 be the subordination functions solving the system 
4’mi,/22(wi(2),W2(z),2;) = 0, z€ C+. 

Then there exists a (small) constant do > 0 such that whenever we have 


\uJi{zo) - wi{zo)\ < do , 

|a;i(zo) - wi(zo)| < do , 

(2.17) 

we also have 

|wi(zo) - uji{zo)\ < 2S'||f(zo)|| , 

|w 2 ( 2 :o) - W 2 (zo)| < 2S'||r(zo)|| • 

(2.18) 


The constant dg > 0 depends on pi and fi 2 , on the interval I as well as on tjm- 


We omit the proof of Corollary 2.6 from Theorem 2.5, since it follows directly from 
Proposition 4.1 in Section 4 below. 

2.3. Applications. We next explain two main applications of the stability estimates ob¬ 
tained in Theorem 2.5. 

2.3.1. Continuity of the free additive convolution. Our first application shows that the free 
additive convolution is a continuous operation when the image is equipped with the topology 
of local uniform convergence of the density in the regular bulk; see (2.23). Bercovici and 
Voiculescu (Proposition 4.13 of [10]) showed that the free additive convolution is continuous 
with respect to weak convergence of measures. More precisely, given two pairs of probability 
measures pLA, and pp on M, the measures pa'^ Pb and pa ffl pp satisfy 

d^ipA ffl PB, Pa ffl pp) < dxipA, Pa) + dx{pB, Pp) , (2-19) 

where du denotes the Levy distance. In particular, weak convergence of pA to pa and weak 
convergence oi pB to pp imply weak convergence of ffl to paW pp. 

Using the Stieltjes transform, we can easily link (2.19) to the systems of equations in (2.5), 
respectively in (2.10). Using integration by parts and the definition of the Stieltjes transform, 
a direct computation reveals that there is a numerical constant C such that 

C / 1 \ 

^ —[l+-jdx{pAWpB,PaWpp) 

C / 1 \ 

< —(1 +-){dx{pA,Pa)+ dx{pB,pp)), 77 = Imz, (2.20) 

Tf \ rjJ 
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for all z € C"*", where we used (2.19) to get the second line. Note that the estimate in (2.20) 
deteriorates as rj approaches the real line. Our next result strengthens (2.20) as follows. We 
consider the measure H “reference” measure (in the sense that it locates the regular 
bulk) while fj-A, IJ'B are arbitrary probability measures and show that the Levy distances 
bound {E + ir/) — {E + ir;)! uniformly in rj, for all E inside the regular bulk 

of Ha ffl /i/3. 

Theorem 2.7. Let Ha o,nd Hi 3 be compactly supported probability measures on K., and assume 
that neither is supported at a single point and that at least one of them is supported at more 
than two points. Let I C be a compact non-empty interval and fix some 0 < t/m < oo. 

Let Ha cLiT'd hb be two arbitrary probability measures on M. 

Then there are constants b > 0 and Z < oo, both depending on the measures Ha o,nd H/ 3 , 
on the interval I as well as on the constant pM, such that whenever 

<ii,{HA,Ha)+ di,{HB,Hp) <b ( 2 . 21 ) 

holds, then 

max |m^^H^g(z) - (z)! < Z(dL(/iA,/ia)+dL(/iB,/i/3)) , (2.22) 

zeSxlo.vM) 

holds, too. 

Note that max^g 52 -(o,r)M) ('^)l < oo by compactness ofl and analyticity of 

in I. Thus the Stieltjes-Perron inversion formula directly implies that {ha H Hb)^^ has a 
density, inside I and that 

< ^(dL(/iA,/ia)+dL(/iB,/i/3)) , (2.23) 

provided that (2.21) holds, where is the density of {ho H 

Remark 2.2. The estimate (2.22) was recently given by Kargin [30] under the assumption 
that (2.14) and (2.15) hold for all z € 5i(0,?7m), i.e., under the assumption that the con¬ 
clusions of our Theorem 2.5 hold. It is quite surprising that one can directly set lm 2 = 0 
in (2.22). As first noted by Kargin, this is due to the regularizing effect of oja, W /3 and to 
the global uniqueness of solutions to (2.5) for arbitrary probability measures. 

2.3.2. Application to random matrix theory. We now turn to an application of Theorem 2.5 
in random matrix theory. Let A = A^^'> and B = be two sequences oi N x N deter¬ 

ministic real diagonal matrices, whose empirical spectral distributions are denoted by ha 
and Hb respectively, i.e., 

1 AT ^ JV 

— ■= , (2.24) 

i—1 i—1 

where A = diag(ai), B = diag(6i). The matrices A and B depend on N, but we omit 
this fact from the notation. Let uja and uib denote the subordination functions associated 
with Ha and hb by Proposition 2.1. 

We assume that there are deterministic probability measures Ha and HP on R, neither 
of them being a single point mass, such that the empirical spectral distributions ha, Hb 
converge weakly to Ha, Hp, bV ^ oo. More precisely, we assume that 

di,{HA,Ha) + di,{HB,Hp)^0, (2.25) 

as N ^ oo. Let uJa, top denote the subordination functions associated with Ha and Hp- 


Let U be an independent N x N Haar distributed unitary matrix (in short Haar unitary) 
and consider the random matrix 

H = := A + UBU* . (2.26) 

We introduce the Green function, Gh, of H and its normalized trace, mn, by setting 

Gh{z) ■■= m_f/(z) := trGi/(z), (2.27) 

z € C+. We refer to z as the spectral parameter and we often write z = E + irj, E £ M., 
77 > 0. Recall the definition of 5i(a, b) in (2.13). We have the following local law for mn- 

Theorem 2.8. Let pa o,nd pp be two compactly supported probability measures on K, and 
assume that neither is only supported at one point and that at least one of them is supported 
at more than two points. Let X C be a compact non-empty interval and fix some 

0 < 77m < 00 ■ Assume that the sequences of matrices A and B in (2.26) are such that their 
empirical eigenvalue distributions pA and ps satisfy (2.25). Fix any small 7 > 0 and set 
Vm :=iV-2/3+T'. 

Then we have the following uniform estimate: For any (small) e > 0 and any (large) D, 

U (2-28) 

^ z^Si{rirn,riyi) ^ 

holds for N > Nq, with some Nq sufficiently large, where we write z = E ip. 

Using standard techniques of random matrix theory, we can translate the estimate (2.28) 
on the Green function into an estimate on the empirical spectral distribution of the matrix iJ. 
Let Ai,..., Aat denote the ordered eigenvalues of H and denote by 

1 ^ 

l^H ■= (2.29) 

i=l 

its empirical spectral distribution. Our result on the rate of convergence of pu is as follows. 


Corollary 2.9. Let X C he a compact non-empty interval. Then, for any Ei < E 2 

in X, we have the following estimate. For any (small) e > 0 and any (large) D we have 


Th{[Ei,E 2)) — PAmB{[Ei, E 2 )) 


> 


fW A 

fV2/3 J 


< N 


-D 


(2.30) 


for N > Nq, with some Nq sufficiently large. 


We omit the proof of Corollary 2.9 from Theorem 2.8, but mention that the normalized 
trace of the Green function and the empirical spectral distribution pn of FI are linked by 


mniz) = tr Gh{z) = 


N 


1 


1 


N ^ X,-z 

i—1 


dpnjx) 
X — z 


zec+. 


Corollary 2.9 then follows from a standard application of the Helffer-Sjostrand functional 
calculus; see e.g.. Section 7.1 of [22] for a similar argument. 

Note that assumption (2.25) does not exclude that the matrix FI has outliers in the 
large N limit. In fact, the model FI = A-\- UBU* shows a rich phenomenology when, say, A 
has a finite number of large spikes; we refer to the recent works in [ 8 , 9, 16, 31]. 


Remark 2.3. Our results in Theorem 2.8 and Corollary 2.9 are stated for U Haar distributed 
on the unitary group U{N). However, they also hold true (with the same proofs) when U 
is Haar distributed on the orthogonal group 0{N). 
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Remark 2.4. In [3], we derive, with a different approach, the estimate (with the notation 
of (2.28)) 

U i |wi?(^) — (^)| > 

^ zeSTirirntriM) '' 

for N > Nq, with some Nq sufficiently large, and with In fact, we obtain 

estimates for individual matrix elements of the resolvent Gh as well. Comparing with (2.28), 
we see that we can choose rj in (2.31) almost as small as at the price of losing a 
factor V^- The stability and perturbation analysis in [3] rely on the optimal results in 
Theorem 2.5 and Theorem 2.7 as well as in Sections 3-5 of the present paper. 

2.4. Organization of the paper. In Section 3, we consider the stability of the system (2.5) 
when at least one of the measures p,i and ^2 is supported at more than two points and we 
give the proof of Theorem 2.5. In Section 4, we consider perturbations of the system (2.5) 
and derive results that will be used in the proof of Theorem 2.8 and also in [3]. In Sec¬ 
tion 5, we prove Theorem 2.7. In Section 6 we consider the random matrix setup and prove 
Theorem 2.8. In the final Section 7, we separately settle the special case when both fii 
and fj ,2 are combinations of two point masses and give the results analogous to Theorem 2.5, 
Theorem 2.7 and Theorem 2.8 for that case. 


3. Stability of the system (2.11) and proof of Theorem 2.5 

In this section, we discuss stability properties of the system (2.11), with p,i, ^2 two 
compactly supported probability measures satisfying the assumptions of Theorem 2.5. 

Lemma 3.1. Let fii, ^2 be two probability measures on K neither of them being supported 
at a single point. Then there is, for any compact set K C C"*", a strictly positive constant 
0 < ( 7 (^ 1 ,/C) < 1 such that the reciprocal Stieltjes transform (see (2.2)) satisfies 


Im2; < (1 - cr(/ri,/C))ImF^j(z), 


yzelC. 


(3.1) 


Similarly, there is 0 < a(p, 2 ,lC) < 1 such that (3.1) holds with p 2 o,nd respectively. 

Assume in addition that pi is supported at more than two points. Then there is, for any 
compact set K, C C+, a strictly positive constant 0 < a(p,i,lC) < 1 such that 


ImT)jj(2:) — Imz 
Imz 


yzGiC, 


(3.2) 


where F(^(z) = d^F^^(z). 

Proof of Lemma 3.1. Assuming by contradiction that inequality (3.1) saturates (with van¬ 
ishing constant a(pi,IC) = 0, for some z € K. C C+), we have Imf)i^( 2 ) = Imz for some z, 
thus (z) = z — a, a G R, i.e., p,i = 5a. This shows (3.1). 

To establish (3.2), we first note that the analytic functions F^. : C’*' —?► C^, j = 1,2, 
have the Nevanliima representations 


(^) = aF,^+^+ f dpF,. (x), 

Jr X — z 


J = l,2, 


z G C+ , (3.3) 


where ap G M and are finite Borel measures on R. Note that the coefficients of z on 

pij F'fj 

the right-hand side are determined by (2.3). From (3.3) we see that 

1 + x^ 


\F’(z)-l\ = 


(x — z)- 


■ (x) 


G C+. 


(3.4) 








10 


as well as 


Im (z) — Im z 
Imz 


1 


■ (x), 


e C+. 


(3.5) 


Hence, assuming by contradiction that inequality (3.2) saturates (with cr(^i,/C) = 0, for 
some z € /C), we must have 


1 + 


X - z 


(x) = 




(3.6) 


for some z € K., implying that Pf^-^ is either a single point mass or Pf^-^ = 0. In the latter 
case, we have (z) = + z and we conclude that pi must be single point measure, but 

this is excluded by assumption. Thus Pf^-^ is a single point mass, i.e., there is a constant 
€ M such that F^^{z) = aF^^ + z + (1 + — z), z £ JC. It follows that pi is a 

convex combination of two point measures yielding a contradiction. This shows (3.2). □ 


3.1. Bounds on the subordination functions. Let pi, p 2 be as above and let wi(z), 
uj 2 (z) be the associated subordination functions. Recall that we rewrite the defining equa¬ 
tions (2.5) for LOi and 1 x 2 in the compact form (wi, 012 , 2 ) = 0 introduced in (2.11). 

We first provide upper bounds on the subordination functions a;i(z), 012 (^). Our proof 
relies on the assumption that pi, p 2 are compactly supported, i.e., that there is a constant 
L < 00 such that 


supp/ii C [—L,L\ , supp ^2 C [—L,L\. (3.7) 

Recall from Theorem 2.5 that we fixed a compact non-empty interval X C Since 

the density is real analytic inside the regular bulk by Proposition 2.3 and since X is 

compact, there exists a constant kq > 0 such that 

0 < Ko < min ( x ) . (3.8) 

x^X 

Fixing a constant 0 < ?7 m < 00 , it further follows that there is a constant M < 00 such that 

max < M. (3.9) 

2e5z(0,?7M) 


Lemma 3.2. Let pi, p 2 be two compactly supported probability measures on M satisfy¬ 
ing (3.7), for some L < 00 , and assume that both are supported at more than one point. Let 
X C be a compact non-empty interval. Then there is a constant K < 00 such that 


max |a;i(z)| < 

zeSx(0,TIM) 


K 


K 


zGox{0,r}M) ^ 


(3.10) 


The constant K depends on the constant rjM o,nd on the interval X as well as on the measures 
pi, p 2 through the constants kq in (3.8) and the constant L in (3.7). 


Proof. We start by noticing that there is a constant ki > 0 such that 


Imm 


/ilffl/i; 


W= / 

JR 


pd{pi ffl P 2 )ix) 


> 


{x — EY -I- ’ 


{x — EY -I- rY 

uniformly hi z = E ip ^ Si{0, pm), where we used (3.8). Thus by subordination we have 

dyii(a) 


min 

zGSx{0,r]M) zGox{0,r}M) 


a — uj 2 {z) 


> Kl , 


(3.11) 
have 

(3.12) 


since TO^iffl^ 2 (z) = -l/^),iffl^ 2 (^) by (2.6). 

On the other hand, pi is supported on the interval [—L, L]; see (3.7). Hence, using (3.12), 
|a; 2 (^:)| must be bounded from above on 5 x(0,?7m)- Interchanging the roles of the indices 1 
and 2, we also get that |a;i(z)| is bounded from above on 5 'x(0,77m)- n 
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Having established upper bounds on the subordination functions, we show that their 
imaginary parts are uniformly bounded from below on the domain 5i(0,7yM)- The proof 
relies on inequality (3.1). 

Lemma 3.3. Let fj,i, fi 2 be two probability measures on K satisfying (3.7), for some L < oo, 
and assume that neither of them is only supported at a single point. Let X C be a 

compact non-empty interval. Then there is a strictly positive constant k > 0 such that 

min Imwi(z)>2fc, min lmuJ2{z)>2k. (3.13) 

zG5z(0,77m) 2eSi(0,7)M) 

Remark 3.1. The constant k in (3.13) depends on the interval I through the constants kq 
in (3.8) and M in (3.9). It further depends on rjM, as well as on a{pi,K 2 ) and cr(/r 2 ,/Ci) 
in (3.1), with 

/Ci = {u e : M = 011 ( 2 :) ,z G 5 i(0,?7m)} , * = 1,2. (3.14) 

Proof of Lemma 3.3. First note that there is ki > 0 such that > ki for all 

z G 5 x(0,?7m); c.f., (3.11). Moreover, there is M < 00 such that < M for all 

z G 5i(0, ryivi); c.f., (3.9). Recall from (2.5) and (2.6) that 

oil (z) + 012 ( 2 ) = 2 ;-yw , 2 G . (3.15) 

Hence, considering the imaginary part, we notice from (3.9) that there is K 2 > 0 such that 

min (Imoii( 2 ) + Imoi 2 ( 2 )) > K 2 . (3.16) 

It remains to show that Im oii and Im 012 are separately bounded from below. To do so we 
invoke (3.1) and assume by contradiction that lmoii( 2 ) < e, for some small 0 < e < K 2 / 2 . 
We must thus have 1111012 ( 2 ) > K 2 / 2 . Since fii is assumed not to be a single point mass. 
Lemma 3.1 assures that 

Imf^^ ( 012 ( 2 )) > , 2 : G 5i(0,7?m) , (3.17) 

1 - a[p.i,K.2) 

with 0 < cr(^i,/C 2 ) < 1, where IC 2 denotes the image of Sz{0,riM) under the map 012 (which 
is necessarily compact by Lemma 3.2). On the other hand, (2.11) implies 

ImF^^ ( 012 ( 2 )) = 1111012 ( 2 ) + Im 011 ( 2 ) — Im 2 , 2 G C'*' . (3.18) 

Since Imoii( 2 ) > Imz, by Proposition 2.1, we get, by comparing (3.18) and (3.17), a con¬ 
tradiction with the assumption that Im 011 ( 2 ) < e, for sufficiently small e. Repeating the 
argument with the roles of the indices I and 2 interchanged, we get (3.13). □ 

3.2. Linear stability of (2.11). Having established lower and upper bounds on the subor¬ 
dination functions oii, 012 , we now turn to the stability of the system (oii, 012 , 2 ) = 0. 

Remember that we call the system linearly S'-stable at ( 011 , 012 ) if ( 011 , 012 ) < S, 

where is defined in (2.12). 

Lemma 3.4. Let /ii, fi 2 be two probability measures on R satisfying (3.7) for some L < 00 . 
Assume that neither of them is a single point mass and that at least one of them is supported 
at more than two points. Let I C be a compact non-empty interval. 

Then, there is a finite constant S such that 

max rfj,^^fj,^{uJi{z),u! 2 {z)) < S, 
z&Sx{0,riM) 


(3.19) 
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and 


max |wi(z)| < 25 , 
^e5r(0.r;M) 


max Iwof-z)! < 25 , 
^eSx(o,w) 


(3.20) 


where uji{z), 012 ( 2 :) are the solutions to ^z) = 0. 


Remark 3.2. Lemma 3.4 is the first instance where we use that at least one of fii and ^2 
is supported at more than two points. For definiteness, we assume that ni is supported 
at more than two points. The constant 5 in (3.19) depends on the interval I through the 
constant kq in (3.8), on the constants L in (3.7), cr(/xi,/C 2 ) and (7(/r2,/C2), as well as on 
5 ( 1 * 1 ,^ 2 ) of (3.2) with IC 2 defined in (3.14). 

Proof of Lemma 3.4- Using (2.12) and Cramer’s rule, F = (oii, 012 , z) equals 


r = 


i-(f;^(o2)-i)(f;Joi)-i)| 


-1 


-F;,(oi(^)) + 1 


-f;^( 02 (^)) + i 


-1 


(3.21) 


As above, we assume for definiteness that is supported at more than two points. We first 
focus on F'^^{ijJ 2 ). Recalling the definition of IC 2 from Remark 3.1 and invoking (3.2), we 
obtain 

- II < (1 ■ (3.22) 

for all 2 ; € 5x(0,?7m), where 0 < cr(/ii,/C 2 ') < 1. Abbreviating a = a(fii,/C 2 ) and using 
^r-i,M 2 (^i(^)^^ 2 (z),z) = 0, we thus have 


lF'(co2(z))-ll<(l-a) 


lma;i(z) 


^ "^lma;2(2;)’ 

Reasoning in the similar way {c.f, (4.9)), we also obtain 


\FMz))-l\< 


Im 012 ( 2 ;) 


z G 5 i(0,77m) ■ 


G C+, 


(3.23) 


(3.24) 


Imo;i(z) ’ 

where the inequality may saturate here since we do not exclude fi 2 being supported at two 
points only. Multiplying (3.23) and (3.24), we get 


max |F;;^(a; 2 (z)) - 1| |U;;^(o;i(z)) - 1| < 1 - cr. 
zGSx{0,'nM) 

Using Lemma 3.2 and Lemma 3.3, we also have from (3.23) and (3.24) that 

K 


max \F'{ujj{z))-l\< 


2e5z(0,r/M) 


4fc 


{*, 7 } = { 1 , 2 }. 


(3.25) 


(3.26) 


Hence, bounding the operator norm by the Hilbert-Schmidt norm in (3.21), we obtain 
by (3.26) and (3.25) that 

1/2 


max J 
2e5z(0,?7M) 


,(a;i(z), 012 ( 2 :)) < 


V2 


1 + 


i-Y 

\AkJ 


=:5, 


(3.27) 


with finite constant 5. This proves (3.19). 

The estimates in (3.20) follow by differentiating the equation ^^ 2 ( 011 ( 2 ), 012 (z), z) = 0 

with respect to z. We get 


-1 

- 1 


- 1 


-1 


AA) 

A{z) 


(3.28) 


From (3.19) we know that 4) is uniformly 5-stable and we get (3.20) by inverting (3.28). □ 
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Remark 3.3. The crucial estimate in the proof above is (3.25). An alternative proof of (3.25) 
under the assumptions that both and 112 are supported at more than two points was 
pointed out by an anonymous referee. From (2.5) we observe that the subordination func¬ 
tion uJi{z) appears, for fixed z G C+, as the fixed point of the map Tz '■ C+ —>■ C+, 

u Tz{u) := (F ^2 {u)-u \ z)- (u) - u + z . (3.29) 

Indeed, assuming that wi G C"*" and that ^ 1 , fi 2 are supported at least at three points 
(so that Ffj,^{z) — z and F^^{z) — z are not Mobius transformations), the fixed point uJi{z) 
is attracting as was shown in [4]. Thus, for any fixed fc > 0, the Schwarz-Pick Theorem 
and (2.5) imply that for any compact subset ^ of {z G C’*' U K : Ima;i(z), Im W 2 (-z) > 2/c} 
there is a constant d{lC) < 1 such that \F{{uj 2 {z) — ~ 1 | ^ any 

z € K,. Thus, under the assumption that and ^2 are both supported at least at three 
points, (3.25) follows from Lemma 3.2 and Lemma 3.3. 

Collecting the results of this section, we obtain the proof of Theorem 2.5. 

Proof of Theorem 2.5. Lemma 3.3 proves (2.14). Lemma 3.4 proves (2.15). □ 


4. Perturbations of the system (2.11) 


In this section, we study perturbations of the system (wi, W 2 , ^) = 0, where p,i, 

[i 2 denote general compactly supported probability measures on R. The main results of 
this section, Proposition 4.1 below, is used repeatedly in the continuity argument to prove 
Theorem 2.8. Yet, as noted in Corollary 2.6, it is of interest itself and it is also used in [3]. 

Proposition 4.1. Fix zq G C"*". Assume that the functions wi, lu 2 , ri, r 2 : C"*" —>■ C satisfy 


Imwi(zo) > 0, Ima; 2 (zo) > 0 and 

^Mi.M2(wi(zo),W2(2:o),2o) =t(zo), (4.1) 

where r(z) := (ri(z),r 2 (z))^. Assume moreover that there is S G [0,1] such that 

|wi(zo) - a;i(zo)| < (5, |w 2 (- 2 o) - W 2 (zo)| < 5, (4.2) 

where coi{z), uj 2 {z) solve the unperturbed system z) = 0. Assume that there 

is a constant S such that is linearly S-stable at {uji{zq),uj 2 {zq)), and assume in addition 
that there are strictly positive constants K and k with k > S and with > SKS such that 

0 < 2fc < Ima;i(zo) < AT, 0 < 2fc < Ima;2(zo) < AT. (4-3) 

Then we have the bounds 

|wi(zo) - a;i(zo)| < 25'||f(zo)||, |w 2 (zo) - a; 2 (zo)| < 2S'||r(zo)|| . (4.4) 

Proof. Combining (4.3) and (4.2) with 6 < k, we get 

Imc5i(zo)>fc, lmuj 2 {zo) > k. (4-5) 


Next, we bound higher derivatives of Fi = , z = 1, 2. We first note that by the Nevanlinna 

representation (3.3) we have 


ImFijui) 

Imw 


= 1-f 


1 + 


a; — w 


■ dPFi (x) , 


G C+, 


z = 1,2. 


(4.6) 


On the other hand, we also have from (3.3) that 

1 -I- a;^ 


|A)'(w)-l|< f dpF,{x), 

Jr \x — <ja| 


w G C+, 


z = 1,2 


(4.7) 
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and analogously for higher derivatives, n > 1 , 

uj € C+, i = 1,2. Thus, combining (4.7), (4.6) and (2.11) we get 


1 + 


(Imtu)’^ 


a; — w 


■ <^PFi (x) , 


(4.8) 


\Fi{oj,iz))-l\< 


lmFi{ujj{z)) — lTaujj{z) Ima;i( 2 ;) — Imz 


hncoj^z) 


linujj(z) 


{z,j} = {l,2}, (4.9) 


eC+, {i,j} = {l,2}. 


(4.10) 


z £ C’*', and similarly, starting from (4.8), 

|p(-)/, . . / Imw,(z)-Imz 

(Ima;,(.))" ’ 

Let ^li{z) := u;i{z) — uJi{z), i = 1,2, and := (ni,n 2 )^- Fixing z = zq and Taylor 
expanding F’i(w 2 (-Zo)) around u! 2 {zo) we get 

1 

nl 


F[(u}2{zo))Vl2{zo) - Oi(zo) - ll2{zo) = ri{z) - V' ^F^'^\uj2{zo))ll2{zoT ■ (4-11) 

Z—/ 7^1 


n>2 


Recalling that ||n(zo)|l/2fc < 6/k < 1 and using (4.10) together with (4.3), we obtain 
from (4.11) the estimate 

\F[{uJ2{zo))1^2{zo) - 111 ( 2 : 0 ) - 112 ( 2 : 0)1 < F(2:o)|| + ^ll^(^o)r , (4.12) 

and the analogous expansion with the roles of the indices 1 and 2 interchanged. We therefore 
obtain from ( 2 . 12 ) and from solving the linearized equation that 

||ll(^o)||<^||r(zo)|| + ^||ll(^o)f. (4.13) 

Thus, we have the dichotomy that either ||ll(zo)|l < 25 ||r|| or 2{KS)~^k'^ < ||ll(zo)|l. Since 
> 5KS by assumption, the second alternative contradicts ||ll(2:o)|| < “25. This proves the 
estimates in (4.4). □ 


In Proposition 4.1 we assumed the apriori bound |wi — Wi| < <5; see (4.2). The next lemma 
shows that we may drop this assumption, for spectral parameters z with sufficiently large 
imaginary part, at the price of assuming effective lower bounds on Imw^. This statement 
will be used as an initial input to start the continuity argument in Section 6. 

Lemma 4.2. Assume there is a (large) rjo > 0 such that for any z £ C’*' with Iraz > rjo the 
analytic functions uji, UJ 2 , ri, ¥2 : C'*' —^ C satisfy 

Imwi(z) — Imz > 2||r(2:)|| , Imw 2 ( 2 :) — Imz > 2||r(z)|| . (4-14) 

and 

^uuU 2 {^i{z),dj 2 {z),z) = r(z), (4.15) 

where r{z) := (ri(z),r 2 (z))^. 

Then there is a constant 770 > 0, with 770 >po> such that 

|wi(z)-a;i(z)| < 2||f(z)|| , |52(2:) - W 2 (z)| < 2||f(z)|| , (4.16) 

on the domain {z £ C"*" : Imz > 770 }, where Wi and 0 J 2 are the subordination functions 
associated with fj,i and fi 2 - The constant rjo depends on the measures fii and p, 2 , and on the 
function r through the constant 770 > 0 . 
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Proof. Recall the Nevanlinna representation ( 3 . 3 ) for and Since /ri and fi2 are 

compactly supported, we have, as Imw ^ oo. 


F^j(w) - w = oi +0(|w| ^), F^2(w) - w = 02 + 0(|w| ^), 

with oi = OiT- and 02 = ap ■ There are thus s^i,S2 : C"*" —> C such that 


<i) 


Ml >1*2 


{UJI{Z),U2{Z),Z) = 


_( ai + si{z) - uji{z) + z \ _ ( ri{z) 


02 + S2{z) — 012 ( 2 ) + 2 


r2{z) 


( 4 . 17 ) 


( 4 . 18 ) 


with 


l,{z) = 0{\^2{z)\-^), ?2(^)=0(Pi(^)|-'), (4.19) 

as lm2; ^ 00. It follows immediately that 011(2) = O(Imz) and 012(2) = O(Imz), as 
lm2 ^ 00. Thus, recalling the definition of in ( 2 . 12 ), we get 

rMi.M 2 (wi(^)i ^2(2)) = 1 + 0(77"^), ( 4 . 20 ) 

as ?7 = Imz 00. In particular, we obtain 

||((D$)-^$)(wi(z),W 2(^),^)|| < ||r^,,^,(wi(2),W2(^))||||4'(o3l(2),W2(2),2)|| 

< 2 ||f( 2 )||, ( 4 . 21 ) 


for Imz sufficiently large. From ( 4 . 8 ) and ( 4 . 17 ), we also get 

01 € C+, 

as Imo; ^ 00. Thus the matrix of second derivatives of $ given by 


Imf 

(Imo;)^ 


D $(011,012) := 1^^^(011,012,2), -^(011,012,2)j 


Fhficzi) 


* = 1 , 2 , 


0 


( 4 . 22 ) 


satisfies ||D^$(oii(2), oi2(2))|| = 0 (Im 2 ;) asimz ^ 00. Hence, choosing 170 > 0 sufficiently 
large, we achieve that 

so :=2||F(z)||||D2$(Si(2),52(z))||<i, 

on the domain {2 G C"*" : Imz > 770}. Thus, by the Newton-Kantorovich theorem 
(see, e.g., Theorem 1 in [ 25 ]), there are for every such z unique oii(z),oi2(z) such that 
4 ‘Mi.M 2 (wi(^),oi 2 ( 2 ),z) = 0 , with 


P^(2)-co.(z)|< ^ ^^° ||f(z)||< 2 ||f(z)||, * = 1 , 2 . ( 4 . 23 ) 

So 

Finally, we note that Imoii(z) = Imoii(z) —Imoii(z)+Imo 3 i > Imz, by ( 4 . 14 ), for Imz > rjQ. 
Similarly, lmoi2(2) > Imz, for Imz > *70. It further follows that r^j_^2(oii(z),o}2(2)) ^ 0 , 
for all z € C'*' with Imz > *70, thus oii(z) and 0)2(2) are analytic on {z € C : Imz > *70} 
since F^^ and F^^ are. Finally, using ( 4 . 17 ) with 01 = wi, 01 = 0)2 respectively, we see that 


lim 

77/^00 ITJ 


lim 

?7/^oo 17 ] 


= 1 . 


Thus, by the uniqueness claim in Proposition 2 . 1 , oii(z),o)2(-z) agree with oii(z),012(2) on 
the domain {z G C+ : Imz > 770}. This proves ( 4 . 16 ) from ( 4 . 23 ). □ 
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5. Proof of Theorem 2.7 


In the setup of Theorem 2.7 we have two pairs of probability measures on M, /Iq,, up and 
MA, MBj where we consider /Iq., M /3 S'® “reference” measures (in the sense that they satisfy 
the assumptions of Theorem 2.7), while ha, Hb are arbitrary. Under the assumptions of 
Theorem 2.7 we can apply Theorem 2.5 with the choices Ha = Hi and Hp = P- 2 - 
Recall from (2.13) the definition of the domain Sx{a, b), a < b. 


Lemma 5.1. Let ha, Hb o,nd Ha, HP probability measures from (2.24) and (2.25) 

satisfying the assumptions of Theorem 2.1. Let u!a,<^b o-nd iOa,eop denote the associated sub¬ 
ordination functions by Proposition 2.1. Let X C be a compact non-empty interval. 

Fix 0 < r]M < oo. 

Then there are a (small) constant bo > 0 and a (large) constant Ki < oo, both depending 
on the measures Ha o-nd Hp, on the interval X and on the constant tjm, such that whenever 


holds, then 


dL{HA,Ha) + dL{HB,Hp) < bo , 


\uja{z) 

|wb( 2 ) 


- U}a{z)\ < Ki 

- ojp{z)\ < Ki 


dL{HA,Ha) 

d-L(HA,Ha) 

{lm.wp{z)Y 


dLiHB,Hp) 

^ (Ima;«(z))2 ’ 

^ dhiHByHp) 

^ (ImWa(z))2 ’ 


(5.1) 


(5.2) 


hold uniformly on 5i(0,77m)- In particular, choosing b <bo sufficiently small and assuming 
that dhiHA, Ha) + dxiHB, Hp) < b, we have, 

max |u;yi(z)| < K, max \ujb{z)\ < K, (5-3) 

zeSx(0,r/M) zeSxiO.VM) 


min ImujA(z)>k, min IinujB(z)>k, (5-4) 

z£Sx(0,'riM) zeSxiO.VM) 

where K and k are the constant from Lemma 3.2 and Lemma 3.3, respectively. 

Remark 5.1. Armed with the conclusions of Theorem 2.5, our proof follows closely the argu¬ 
ments of [30]. We further remark that the main argument in the proof of Lemma 5.1 is differ¬ 
ent from the ones given in Section 4: it crucially relies on the global uniqueness of solutions 
on the upper half plane for both systems, [uia,cjp, z) = 0 and {uja,cjb, z) = 0 , 

asserted by Proposition 2.1. 


Proof of Lemma 5.1. We first write the system (wq,(z), w^(z), z) = 0 as 


<i> 




(uJa(z),ujp(z),z) = r{z), 


G C+, 


with 


/ \ - ( rAiz) 

V '^b{z) 


(5.5) 


._ f FnA{oJp{z))-F^^{u}p{z)) 

\ F^^{uJa{z)) — Ffj_p(uJct{z)) 

From Lemma 3.3, we know that the imaginary parts of the subordination functions uja,ujp 
are uniformly bounded from below on iSx(0, tjm)- Next, integration by parts reveals that for 
any probability measures Hi a-nd H 2 , 


|m^,(z) - mf,fz)\ < c 


dL(A^l,M 2 ) 

Imz 


1 + 


Imz 


with some numerical constant c; see, e.g., [31]. Thus, 


\Ff,^{ujp{z)) - Ff,^{ujp{z))\ = 


I™A»a(^/ 3(^))) - (^/9(^))l 

\m^^{ujp{z))mf,^ {cop{z))\ 


< C 


ZGC+, 

dx{HA,Ha) 

{lm.ujp{z))'^ 


(5.6) 


(5.7) 
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with a new constant C that depends on the lower bound of Imm^^(w/ 3 ( 2 )) = 
which is strictly positive on 5 i(0,?7m); c./., (3.12). Here we used 


lTam^^{u:p{z)) > lTtvm^^{LOp{z)) - \lmm^j^{uip{z)) - Imm^„( 01 / 3 ( 2 ))| > ( 01 ^( 2 )), 

as follows from (5.6) for small enough ^^a) < 5o- Repeating the argument with the 

roles of A and B interchanged, we arrive at 


\rA{z)\ < C 


dL(MA,Ma) 


\rB{z)\<C 




zGSi{0,r]M), (5.8) 


{luiUJp[z)Y' I'-OV";! (luiOJaiz))'^ ' 

for some constant C. Recalling the definition of T in (2.12), we get for sufficiently small 60 , 

(oiq,, 01/3) ^ 2 r^^ (oict, 01/3) ^ 2 * 5 , ( 5 . 9 ) 

where S is from Lemma 3.4, and where we also use Lemma 3.3 and the assumption di^{fiA, Ma) + 
^hifJ-B, fJ'p) "£ bg. The Newton-Kantorovich theorem then implies {c.f., the proof of Lemma 4.2 
for a similar application) that there are 01 ^( 2 ), oj_b( 2 ) satisfying 


(wa( 2 ), UJBiz), 2 ) = 0 , 


2 € 5i(0,?7m) , 


(5.10) 


and 


\uJa{z)-QA{z)\<2\\r{z)\\, 101 / 3 ( 2 )-o;b(^)| < 2||r(2)|| , (5.11) 

2 G 5x(0,77m). Invoking (5.8), (5.11) and Lemma 3.3, we see that lja{z) G C+ and 015 ( 2 ) G 
C+, for any 2 G 5i(0, 7 ?m) if 5o is sufficiently small. Yet, by the global uniqueness of solutions 
asserted in Proposition 2.1, we must have ujAiz) = 01 ^( 2 ), ljb(z) = 015 ( 2 ), 2 G C+. Together 
with (5.11) and (5.8) this implies (5.2) and concludes the proof. Then, choosing bo sufficiently 
small, (5.3) and (5.4) are direct consequences of (5.2), Lemma 3.2 and Lemma 3.3. □ 


With the aid of Lemma 3.4, we prove the stability of the system (oia, 015 , 2 ) = 0. 

Corollary 5.2. Under the assumptions of Lemma 5.1, there is a (small) constant bi > 0, 
depending on the measures fia cind pLp, on the interval I and on the constant rjM, such that 

dL(MA,/ia)+dL(AiB,/i/3) < 5i (5-12) 

implies 


max Tfj,^^fj,^{ujA{z),uJB{z)) <2S 

zeSz(0,TlM) 


(5.13) 


and 


max |oi(i( 2 )| < 4S', 
ze<Si(o,7?M) 


max 1015 ( 2)1 < 45', 
ze5x(0.r)M) 


(5.14) 


where iiJA{z), ujb{z) satisfy ^^^^^g{iOA{z),uJB{z),z) = 0 and S is the constant in Lemma 3.4- 
Proof. Let T = ( 01 ^( 2 ), 015 ( 2 )). Analogously to (3.21), we have 


r = 


1 




-1 


-F(Juja{z)) + 1 


-Ff^iwB{z)) + l 


-1 


(5.15) 


Using the bounds (5.3) and (5.4) for sufficiently small 61 , we follow, mutatis mutandis, the 
proof of Lemma 3.4 to get (5.13). The estimates in (5.14) then follow as in Lemma 3.4. □ 


We are now ready to complete the proof of Theorem 2.7. 
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Proof of Theorem 2 . 1 . Recall that (013(2)), 2 € C+. We first note that 


mfj,^{uJB{z)) - m^„(wB(^))| < C 


dL(A^A,/la) 
Im 01^(2) 


1 + 


Im uJ/siz) 


for some numerical constant C; c.f., (5.6). Thus using (5.4) we get 

\mf_c^{ujBiz)) - mfj_^{ujBiz))\ < K2k~‘^dBifJ-A, h'a) , Z e 5 i(0,77m) , 

for some numerical constant K 2 . Choosing b as in Lemma 5.1 and assuming that h'a)+ 

dhifJ-A, h'is) ^ '"^6 get from (5.2) that 

|w^„(wb( 2 )) - 771^^(01,3(2)) I < Rrifc“^((dL(MA,/ia) + d^iHA, h'l^)) , Z G 5 x( 0 , 77 m) . 

Setting Z := Kik~‘^ + K 2 k~^ we thus obtain (2.22). □ 

Remark 5.2. Note that under the assumptions of Theorem 2.7, we have, for d^i^-A, h'a) + 
dL(/7A,/i/3) < b, the bounds 

ki/ 2 < < 1/fc, (5.16) 

uniformly on iSx(0, 77 m) with 7ti > 0 from (3.12) and k > 0 from (5.4). 

6. Proof of Theorem 2.8 

Before we immerse into the details of the proof of Theorem 2.8, we outline how The¬ 
orem 2.5 and the local stability results of Section 4 in combination with concentration 
estimates for the unitary groups lead to the local law in (2.28). 

6.1. Outline of proof. We briefly outline of our proof when U is Haar distributed on U (TV). 
Since we are interested in the tracial quantity mn oi H = A + UBU*, we may replace H 
by the matrix 


H := VAV* + UBU* 


( 6 . 1 ) 


where V is another Haar unitary independent from U. By cyclicity of the trace we have mu = 
and we study below. We emphasize that this replacement is a convenient techni¬ 
cality which is not essential to our proof. 

Using the shorthand 


G C+ 


A := VAV* , B := UBU* , 
we introduce the Green functions 

Gj(z):=(A-z)-\ Gg(z) := (B - z)-\ 

For a given N x N matrix Q, we introduce the function 

fqiz) :=trQG^{z), 2 G C+ 

where G^ = {H — z)~^ is the Green function of H. We define the approximate subordination 
functions, uj\ and W3, by setting 

E/a(^) 


( 6 . 2 ) 

( 6 . 3 ) 

( 6 . 4 ) 


uj\{z) := 2- 


w|j(2) := 2 - 




G C+, 


( 6 . 5 ) 


M,mpj{z) ’ M,mjj{z) ’ 

where the expectation E is with respect to both Haar unitaries U and V. From the identity 
{H — z)Gjj{z) = 1 , 2 G C+, we then obtain the relation 

1 


Uj'Xiz) + UJ%{z) - 2 = - 


Em^{z) ’ 


G C+, 


( 6 . 6 ) 
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reminiscent to (c./., (2.5)-(2.6)) 

uja{z) + u}b{z) - z = - 


:ec+. 


mASBiz) ’ 

For the proof of Theorem 2.8, we decompose 

mij{z) - m.AmB[z) = {nijjiz) - Em^{z)) + {Em^{z) - mAmB{z)) 


(6.7) 


where we abbreviate To control the fluctuation part, m^{z) — E'm^{z), 

we rely on the Gromov-Milman concentration inequality [26] for the unitary group; see (6.22) 
below. To control the deterministic part, we first note that, by ( 6 . 6 ) and tTiAmB^z) = 
iTiAiujBi{z)), bounding |Em^(z) — 'mAmB{z)\ amounts to bounding |a;^( 2 ) — ijJa{z)\ and 
\uJb{z) — ujb{z)\. We then show that uj'^{z) and ujg{z) are both in the upper-half plane and 
satisfy 


=r{z), z e Sx{r]^,r]M), ( 6 . 8 ) 

for some small error r{z) G C+, i.e., we consider ( 6 . 8 ) as a perturbation of the sys¬ 
tem <l>fj_^^f^^(oJA(z),uJB(z),z) = 0; c./., (2.10). The formal derivation of ( 6 . 8 ) goes back 
to Pastur and Vasilchuk [33]. Using Proposition 4.1 (with rough a priori estimates on 
\uj\{z) — uja{z)\ and \ujg{z) — ujb{z)\ obtained from the continuity argument below) and 
stability results of Theorem 2.5 and of Section 5, we then bound \uj‘^{z) — oja{z)\ and 
\u!g{z) — ijJb{z)\ in terms of r{z). 

In sum, for fixed z G C^, our proof includes two parts: (i) estimation of the error r(z) 
in ( 6 . 8 ) and (ii) concentration for m^{z) around Em^( 2 :). Both parts rely on the estimates 

Em^{z), uj'xiz), ujb{z) ^ 1, lmuj‘x{z),lmuj%{z)>l, z G Sx{r]„„r]M) ■ (6.9) 


Note that the quantities in (6.9) are obtained from the Green function of H by averaging 
with respect to the Haar measure. Similar bounds for niAmB, and uib were obtained in 
Section 5. These latter quantities are defined directly from fj,A and ^b via Proposition 2.1 
To establish (6.9), we use a similar continuity argument as was used for Wigner matrices 
in [24]: For Imz = rju sufhciently large, the estimates in (6.9) directly follow from definitions. 
For z = E + ir], with E G I fixed, we decrease rj = rju down to rj = r/m in steps of size 
0{N~^), where, at each step, we invoke parts (i) and {ii). However, a direct application of 
the Gromov-Milman concentration inequality for part (i) does not allow to push 77 below 
the mesoscopic scale 77 = Indeed, the Gromov-Milman inequality is effective if 

jN = 0 ( 1 ), where £ is the Lipschitz constant oim^{z) with respect to the Haar unitary V. 
It is roughly bounded by y^tr jG^(z)j'*/iV, which in turn is trivially bounded by 
giving the 77 > 7 > 0, threshold. However, in reality, the random quantity 

Y^tr \Gjj{z)\'^/N is typically of order 1/^jNrf as follows by combining the deterministic 
estimate tr \G^{z)\^ < ri~^lmmjj{z) with a probabilistic order one bound for Imm^(z). 

Our key novelty here is to capitalize on this latter information. We introduce a smooth 
cutoff that regularizes mjj{z) and then apply the Gromov-Milman inequality for this regu¬ 
larized quantity. With the bound l/y^Nrj^ for the Lipschitz constant, we get concentration 
estimates down to scales 77 > 7 > 0. 


6.1.1. Notation. The following notation for high-probability estimates is suited for our pur¬ 
poses. A slightly different form was first used in [21]. 

Definition 6.1. Let 

X = {X^^\v) : N gN,vGV^^^), 


r = (yW(7;) : Af G N, u e VW) (6.10) 
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be two families of nonnegative random variables where is a possibly N-dependent pa¬ 

rameter set. We say that Y stochastically dominates X, uniformly in v, if for all (small) 
e > 0 and (large) D > 0, 



for sufficiently large N > NQ{e,D). If Y stochastically dominates X, uniformly in v, we 
write X <Y. If we wish to indicate the set explicitly, we write that X{v) -< Y{v) for 
all V G VW. 

6.2. Localized Gromov Milman concentration estimate. In this subsection, we de¬ 
rive concentration bounds for some key tracial quantities. They are tailored for the continuity 
argument of Subsection 6.3 used to complete the proof of Theorem 2.8. The argument works 
with U, V independent and both Haar distributed on U{N) or on 0{N). Below, E denotes 
the expectation with respect Haar measure. 

In the rest of this section, we let I C denote the compact non-empty subset fixed 

in Theorem 2.8. Also recall from Theorem 2.8 that we set rjm = 7 > 0. Below we 

choose the constant tjm ~ 1 to be sufficiently large at first, but from the proof it will be clear 
that we can eventually choose 0 < t^m < 00 arbitrary. Recall from (6.4) the notation /q, 
where Q is an arbitrary N x N matrix. 

Proposition 6.2. Let Q be a given N x N deterministic matrix with ||Q|1 < 1. Fix E G I 
and T] G [ 77 m, Pm] ■ Then 

lmmpj{E -{-ip) ^1, ^P^[p,Pm], (6.12) 

implies the concentration bound 

\fvQv* {E -{- ip) - EfvQv- {E + ip)\ F . ■ (6.13) 

yj N‘‘p'^ 

The same concentration holds with VQV* replaced by UQU*. 

Proof. For fixed E G I, we consider z = E -{- ip G C“'" as a varying spectral parameter and 
use z = E -{-ip for the specific choice in the lemma. By the definition of /(.) and cyclicity of 
the trace, we have 

fvQv-iz) =tTVQV*{VAV* + UBU* - z)~^ =tTQ{A + V*UBU*V - z)~\ (6.14) 

where tr (•) stands for the normalized trace. For simplicity, we denote 

W:=V*U, H:= A + WBW*, Gu{z) := {% - z)-^. (6.15) 

Observe that W is Haar distributed on U{N), respectively 0{N), too. By cyclicity of the 
trace we have tTG-uiz) = tr Gh{z) = mniz). According to (6.14) and (6.15), we may regard 
in the sequel fvQV* as a function of the Haar unitary matrix W by writing 

h{z) = hw{z) := fvQV-{z) • 

For any fixed (small) e > 0, let £ be a smooth cutoff supported on [0, 2N^], with x(x) = 1, 
X G [0, A^^], and with bounded derivatives. Since mniz) = tr Gniz), we can regard mniz) 
as a function of W and write 


X{z) = Xwiz) := xiimmHiz)). 


( 6 . 16 ) 
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We then introduce a regularization, hw, of hw by setting 

1- log2 'll 

h{z) = hw{z) ■= hw{z) xw{E+ i2'^vi). (6.17) 

n—0 

We will often drop the W subscript from the notations hwiz), hw{z) and xw{z) but 
remember that these are random variables depending on the Haar unitary W. 

We will use assumption (6.12) at dyadic points, i.e., that 

luirriHiE + i2^^) ^ 1, 0 < Z < [—log2 iy] , (6.18) 


(recall that mniz) = m^{z) so we may drop the tilde in the subscript of to). Hence, 
by (6.16) and (6.18) we see that, for arbitrary large H > 0, 


1- log2 'll 

n x{E+ 12’^ ifj = 1, i.e., h{z) = h{z), (6.19) 

1=0 

with probability larger than 1 — N~^, for N sufficiently large (depending on e and D). 
Taking the trivial bound IIQII/iy for Zz,(z) and for h('z) into account, we also have 

Eh{z) - E h{z) = O {N-^+^) . (6.20) 


To prove (6.13), it therefore suffices to establish the concentration estimate 

1 


h{T) — Eh{z) 


-< 




( 6 . 21 ) 


for the regularized quantity h(z). 

To verify (6.21), we use the Gromov-Milman concentration inequality [26] (see Theo¬ 
rem 4.4.27 in [2] for similar applications) which states the following. Let Ai{N) = SO(iV) 
or SU(A^) endowed with the Riemann metric ||ds ||2 inherited from Miy[{C) (equipped with 
the Hilbert-Schmidt norm). If g : ||ds|| 2 ) —R is an £-Lipschitz function satisfying 

Eg = 0, then 

IP(l5l > ^) < , V5>0, (6.22) 


with some numerical constant c > 0 (depending only on the symmetry type and not on N). 
Here P and E are with respect Haar measure on Ai{N). 

In order to apply (6.22) to the function W hwi^) = h{z), we need to control its 
Lipschitz constant. To that end, we define the event 

^{rj) = ■= {irnmniE + i2"?7) < 2iV® : Vn G No} . (6.23) 

To bound the Lipschitz constant, we need to bound quantities of the form tr \ G'n{'z)\^ re¬ 
stricted to the event Let (Ai('H)) denote the eigenvalues of TL and introduce 

In:=[E-2''^rj,E + 2^rj]nl, := |{i : A,;(H) G I„}|, n G No . 


Since and 'H are unitarily equivalent, their empirical eigenvalue distributions are the 
same, hh', c.f., (2.29). Using the definition of the Stieltjes transform we have, for all n G No, 
the estimate 

Nr, = N f dgH<lN- 2"+iiy ’ 2”gImTOff(2). 

Jx„ 2 7i5-2"5) (a; - Ey + 

Thus we have on the event H(7y) that 

7V„ < 2'^N^+^rj, 


Vn G No 


(6.24) 
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By the spectral theorem, we can bound 


N 


tr|G«(z)| 


1 


N^^\x,{n)-E\ + fi- 

Then we observe (with the convention I_i = 0) that 

1 1 ^ 1 


(6.25) 


N , ^ AT oo 

N ^ \K{n) ^ In \I„-l) 

fclogA/] 

= V E E 


lx,(n)-Ei + y 


l\(n)-Ei + ^ 


n=0 AigI„\I„ 

where we used ||'H|| < C to truncate the sum over n at [clogA^]. We then bound 


fclogWl 

E E 


r-^iogivi 

E 5;;|Siv'iogiv, 


1 


n=0 Aiel„\l, 

where we used (6.24), i.e., with (6.25) we arrive at 

l(f2(r7))tr|G«(z)| <iV^logiV. 
Using the spectral decomposition of 'H we see that 


1 " 


1 


= ^y- 


N 




N ^^\\{n)-E\^ + rp Nri^^\X,in)-E\^ + rp 


Im rriH^z) 
rf 


(6.26) 


, (6.27) 


where we also used that trG'H(I) = trGi/(S") = tojy(z). Thus, we bound 

I(f 2 (r 7 ))tr |G-h(z)|'= < l(n(r}))rr’'+^lmmH{z) < , Vfc > 2. (6.28) 

Having established (6.26) and (6.28), we proceed to estimate the Lipschitz constant of h{z) 
as a function of W. Let su(iV) and so(iV) denote the (fundamental representations in Mm{C) 
of the) Lie algebras of SU{N) and SO{N) respectively. Let m stand for either su{N) or 
so(A^). Note that X S m satisfies X* = —X. Since SU{N) and SO{N) are matrix groups 
the Lie bracket of su(A^) and so{N) respectively is given by the commutator in the matrix 
algebras. For fixed X G Mn{C), we let adx ■ XIn{C) —>■ Mn{C), Y >->• a,dx{Y) := 
XY — YX . For G m and t G M, we may write (WBW*) = (e*^W)B{e*^W)* , where 
we used that X* = —X. Further note that 

^y^TBW*) = e^^'^^eidxiWBW*) . (6.29) 

For X G m with ||X ||2 = 1, we let 

Gn{z, tX) := (h + {WBW*) - 2) teR, 

and denote accordingly 

mH{z,tX) := tTGn{z,tX) , x{z,tX) := x(Imm//(z,tX)), 

1- log2 »?1 

h{z,tX) :=tTQGn{z,tX) x(F; + i 2 '? 7 , tX), 

1=0 


with x( 2 :, 0 ) = x{z), h{z,0) = h{0), etc. 
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Evaluating the derivative of h{'z,tX) with respect to t at t = 0 we get 
f) ^ 1 - log 2 

-h{z,tX)\^^ = -tT(QGH{z)i^dx{WBW*)Gniz)) x{E + ^2% 

* ° 1=0 
/ F- log2 V^ \- log2 V^ 

-tv{QGn{z))[ Y. f n x{E + y2%\-cj>{E + i2^^) 


i=o 


1=0 


X Imtr(G«(£; + i2^?7)adx(VEBlE*)G«(£; + i2^?7)H , (6.30) 


where we used (6.29) and where we introduced 4>{z) := x^(Inim//(z)), with ^ the derivative 
of X. Recalling (6.16) and the definition of the cutoff x, we note the bounds 
F- log2 ^1 F- log2 

n x(i^ + i2'?y) < 1, [ n x{E + i2'fy)] • cj^{E + i 2 ^i)) = 0(log N ). (6.31) 

1=0 i=o 

On the event the complementary event to ^(iy), we further have the identities 

F- log2 ^1 F- log2 ^1 F- log2 


n x(i? + i2'??) = 0, ^ 

i=o 


x{EEV2}?i) ■(j){E + iVrj) = Q. (6.32) 


1=0 


;=o 


It thus suffices to bound (6.30) on the event n(? 7 ). We bound the first term on the right side 
of (6.30) as 


F- log2 »?1 


V{Sl{^))\iT(QGH{z)a.Ax{WBW*)Gu{z)) J] x{EEi2% 


1=0 


F- log2 ^1 


<l(r!(57))^||adx(lEi?IE*)||2||G«(£)gG«(z)||2 x{E + V2}ri), (6.33) 


1=0 


where we used cyclicity of the trace and Cauchy-Schwarz inequality. Next, note that 
||adx(IERIE *)||2 < 2 ||R|| ||X ||2 < 2||i3||, where we used the definition of adx, ||IE|| < 1 and 
iix ||2 = 1. Similarly, we have \\Gn{z)QGH{z )\\2 < \\Q\\\\Gh{z)G’!^{z)\\ 2 . Thus from (6.33), 

F- iog2 v^ 

l{nm\tT{QGn{z)adxiWBW*)Gniz)) [] xiE + i2^r}) 


1=0 


<2||R||||g||l(fl(iy)) 


tr\Gniz )\ 
N 


4 \ 5 r-log2?l 


n x{E + i 2 '? 7 ) 


i =0 


< 


Nr}^ 


(6.34) 


where we used (6.28) with A: = 4 in the last step. 
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To handle the second term on the right side of (6.30), we use (6.31) and (6.26) to get 
1- log2 ^1 1- log2 VA 

l(0(r))) tr(QG«(z)) ^ J] x{E+ ^2% ■ cj>{E+ 


3=0 


1=0 

1^3 


X Imtr(G-H(£; + i2^77)adx(VF.BlT*) G-h(G + 12^'?))) 

1- log2 'tl 1- log2 

<l(0(57))||Q||tr|G«(z)| ^ J] xiE + i 2 ^rj) + i 2 %\ 


1=0 

1^3 


X 2 ||i?||||g|| 


3=0 

1-^3 

tr|G«(£; + i2J'77)|4\^ 
N 


< 




\Nrp J 

Combining (6.35) and (6.34) we obtain, for any X € su(N) or so(N) with ||X ||2 = 1, that 

•Ar4e 


p(Ztx) 


t=0 


< 


Nr}- 


(6.35) 
1 , that 

(6.36) 


i.e., the Lipschitz constant of h{z) as a function of W is bounded by G 
constant G depending only on ||i?|| and ||(5||. Thus, taking 

'N^e\i 


1/2 


for some 


g = h{z) -Eh{z) , 


C = C 


Nfj- 


6 = N^yy/N^r}^ 


in (6.22), and choosing £ > 0 sufficiently small, we get (6.21). Together with (6.19) and (6.20) 
this implies (6.13). □ 


6.3. Continuity argument. In this subsection, we often omit z € C'*' from the notation. 
Let U and V be independent and both Haar distributed on either U(N) or 0{N). Recalling 
the notation in Section 6.1, we set 

Aa^z) := -(K[mHiz)])G^iz) - (E[/g(z)])G^(z)G5(z), 

AB{z)-.= -mmHiz)])Gj^iz)-mfyz)])Gyz)Gj^{z), z G C+, (6.37) 

where we introduced EX := X —EX, for any random variables X. Using the left-invariance 
of Haar measure, one derives the identities 

E[Gjj 0 IGjj] = E[AG^ 0 Gjj], E[Gjj 0 BG^] = E[BG^ 0 G ^]; 

see Theorem 7 in [33] or Appendix A of [31] for proofs. Taking the partial trace for the first 
component of the tensor products, we get 

EG^iz) = EGy^Biz)) + Syz) , Syz) := ^^-^E[G^(w^(z))(A - z)AAiz)] , 

(6.38) 

where u!%{z) is defined in (6.5), we used ( 6 . 6 ), and where we implicitly assumed that 
Ima;^(z) > 0. This last assumption will be verihed along the continuity argument. Then, 
we set 

_ ^rSyz) _ 

trGyuj%{z)){trGyyyz)) + trS'x{z)) ’ 


ryz) : 


(6.39) 
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and define SU^) and rU^) in the same way by swapping the roles of A and B. Us¬ 
ing (6.38), (6.6), we eventually obtain, under the assumption that > 0 and 

l\nu:%{z) > 0, 

^a‘a.mb(wa( 2 :):W|( 2 :),z) = r“(z), z e C+, (6.40) 

withr=( 2 ) = {r\{z),r%{z)y. 

Lemma 6.3. Fix E G T and any rj G [rj-an Vm] ■ Set the notation z = E + irj and z = E + iri. 
Suppose that 


\uj‘X{z) - uja{z)\ + \ujb{z) - ujb{z)\ < N , yp = lmz G [p,r]M] ■ (6.41) 

Moreover, assume that for the event 

5 ( 77 ) = S£;(r)) := ||TOi/( 2 :) - to^Hb( 2 :)| < : z = E + \r], Vry £ [?),? 7 m]| 


we have 

P(S(77)) > 1 - 7V-^(1 + N^pm - V)) , 
for any D > 0 if N > Ni{D). 

Then, for any e > 0, the estimates 


yAZ)\ + \r%m<y^,. 

\uj'X{z) - ujAiz)\ + \ujb{z) - u ;_ b ( z )| < , 

iV" 

\EmH{z) - mASBiz)] < , 

hold for any N > tV 2 (e). Moreover, for any e,D > 0, the event 

r '1 

0 ( 77 ) = 0e(77) := EEiv) n yruHiz) - mASiB{z)\ > 


(6.42) 

(6.43) 

(6.44) 

(6.45) 

(6.46) 


satisfies 

P(0(i7)) < N-^ , (6.47) 

if N > N 3 {e,D). The threshold functions Ni, N 2 , N 3 depend only on p,a,TP 7 Sie speed of 
convergence in (2.25) and they are uniform infjG [ 77 m, 77 m] and E Gl. 


We postpone the proof of Lemma 6.3 and prove Theorem 2.8 first. 


Proof of Theorem 2.8. We start with observing that it is sufficient to prove a version of (2.28) 
where the real part of the spectral parameter E is fixed. This version asserts that there is a 
large (iV-independent) 77 M, to be fixed below, such that for any (small) e > 0 and (large) D, 
and any fixed E G T, 

U {hH{z)-m^^si^Jz)\ > P/ 2 I) ( 6 - 48 ) 

^ zeSEiVn^.VM) ^ ^ 

holds for N > Nq, i.e., the set 5i(77in,77M) in (2.28) is replaced with 5 b( 77 m, 77 m) := {i?-l-i77 : 
77 G [r]m, 77m]}- The threshold Nq depends on e, D, T and on the speed of convergence 

in (2.25). 

Indeed, by introducing the discretized lattice version 

5i(a, b) := Sx{a, h) fl N~^{1 x iZ} 
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of the spectral domain Sx{a, b) {c.f., (2.13)) and by taking a union bound, we see that (6.48) 
implies 




|TOi/(z)-TO^^ffl^ 3 (z)| > 




7V(Im z)3/2 


< 


C 


ND-5 ■ 


(6.49) 


Thanks to the Lipschitz continuity of the Stieltjes transforms mH{z) and with 

Lipschitz constant rj~‘^ = (Imz)“^ < for any Imz > rj^, we see that (2.28) follows 
from (6.49) after a small adjustment of e and D that were anyway arbitrary. 

From now on we fix i? G I and our goal is to prove (6.48). We will use Lemma 6.3. In 
the first step we verify that the assumptions of this lemma hold for rj = r]M, i-e., that (6.41) 
and (6.42) hold for 2 = + iTyM- In the second step, we successively use Lemma 6.3 to 

reduce rj steps by steps of size N~^ until we have verified (6.41)-(6.42) down to rj = r]^. 
Then (6.48) will follow from a final application of Lemma 6.3 combined with discretization 
argument similar to the one above, but this time the 77 variable instead of the E variable. 

Step 1. Initial bound. First we note that since pA and fiB are compactly supported, ||iL|| 
is deterministically bounded, we thus have !mmH{E + iriM) < (? 7 m)~^ ^ 1 assuming t^m > 1 - 
Following the main argument in the proof of Proposition 6.2, we have the concentration 
inequality 

\fvQV‘{E + ipu) - EfvQV‘{E + i77M)| ^ , 3 , (6.50) 

uniformly for any deterministic Q with ||(5|| < 1. The analog concentration holds with V 
replaced by U. Using (6.50) with Q = I (I the identity matrix), we have | W.[mH{E+ir]M)]\ -< 
N~^. Hence, it suffices to show that 

1 


lEmniE + iriM) - rnABBiE + iriM)\ -< 


N ■ 


(6.51) 


Recalling the definitions of and ujg in (6.5), we have, with 2 = E + ir/M, the expansion 

0 ( 2 -^), 


EtrAG'77(2) trH2 ^+EtT:A{A + B)z 

UJ — z ^ — z 


EtrGg( 2 ) 

as 77 m ^ 00 . Thus using the assumption tr H = 0 we get 


Ima;;^(£' + i77M) - Im77M = 


tr H^77 m + E tr ABtjm 


+ 0 


Vm 


\E + i77MP 

as 77 m 00 . Next, since V and U are independent, we have 

Eti-V AV*UBU* = trE[VAV*]E[UBU*] = trHtrH = 0, 
since tr H = tr H = 0 by assumption. Thus 


Ima;^(F; + i77M) - Im77M = 


tr A^ij 


M 


o 


1 

Vm 


(6.52) 


\E + i77MP 

as 77 m 00 . Since tr > 0, we achieve by choosing rju sufficiently large (but independent 

of N) that 

1 trA^rju 


Imuj'XiE + iriM) -ImriM > 


(6.53) 


2 |U + i77MP ’ 

and the analogue estimate holds with A replaced by B. In particular, we have, for such 77 M, 
Ima;)j(F; + i77M) > 1, ImwKU + i77M) > 1, (6.54) 

and uj'xiE + i77M) -- 1, 0 }%{E + i77M) -- 1. 
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To show (6.51), we apply Lemma 4.2 to the system (6.40). Having established (6.53), it 
suffices to show that 

IrliE + ir]M)\ < ^ , IreiE + ir]M)\ < , (6.55) 


since then we have, for N sufficiently large and 77 M as above that, for any fixed e G [0,1), 


^ ^ 1 

N - 2|S + iryMP 


< Imuj'XiE + iriM) 


Im?7M , 


(6.56) 


and similarly with B replacing A. In particular, combining (6.55) and (6.56), we see that 
assumption (4.14) of Lemma 4.2 (with the choice r = r°) is satisfied for N sufficiently large 
(with high probability). Consequently, we see that (6.41) (even with instead of N~'^ 

in the latter) hold ioi z = E + irjM- Finally, the equations 


EniHiz) 


1 

Z - ^%{z) ’ 




1 

Z - UJAiz) - UJb{z) ’ 


(6.57) 


together with the concentration estimate (6.50) yield (6.42). 

It remains to justify (6.55). Since irjM m 7 f( 7 ?+ i? 7 M) = 0(^m )’ have EmH{E + iriM) ~ 
1 . In addition, from (6.54) it follows that mA{uJB{E+ ir]M)) ^ 1 and mB{u}'^{E + iijm)) ^ I. 
Recalling the definition in (6.39), we see that (6.55) is equivalent to 

S'xiE + irjM)\ ^ ^ , |tr^B(£' + i??M)H ^ ■ (6.58) 

By the definitions of <5^, in (6.38), and A^, in (6.37), it is easy to obtain (6.58) by 
using (6.50) and Cauchy-Schwarz. This completes Step 1, i.e., the verification of (6.41)- 
(6.42) for 77 = r]M- 

Step 2. Induction. Recall that loa, ojb and toabb (see Lemma 5.1) are uniformly bounded 
and uj^{z), ijjg{z), mH{z), loa(z), ujb(z) and myi^B(z) are Lipschitz continuous with a 
Lipschitz constant bounded by (Im^)"^ < N^, for any Im^; > Applying Lemma 6.3 to 
conclude (6.44) with the choice e = 7 /IO, we see that if (6.41) and (6.42) hold for some rj, 
then (6.41) also holds for rj replaced with rj — N~^ as long as rj > t]^. Moreover, by the 
Lipschitz continuity of ttih and mAaB, notice that 


E{rj- N ®) D ^{p)\e{p). 


(6.59) 


Thus, if (6.42) holds for some rj, then (6.59) and (6.47) imply that (6.42) also holds for 77 
replaced with rj — N~^. Using Step 1 as an initial input with the choice rj = tjm, and 
applying the above induction argument O(N^) times by reducing 77 with stepsize N~^, we 
see that (6.41) and (6.42) hold for all rjk G [ 77 ^, rju] of the form rjk = rju — k-N~^ with some 
integer k. Applying Lemma 6.3 once more for these rjk, but now with an arbitrary e > 0, 
we conclude from (6.42) and (6.47) that 


\mH{E + i?7fc) 


mAmB[E + i?7fc)| ^ ^ , 


fc = 0,1,..., fco , 


(6.60) 


where fco is the largest integer with jjkg > rj^. The uniformity of (6.60) in k follows from the 
fact that the threshold functions Nj in Lemma 6.3 are independent of 77 . Clearly fco = 0{N^), 
so taking a union bound of (6.60), compensating the combinatorial factor CN^ by replacing 
D by D—5, and slightly adjusting e to extend the control from the set {z = E+\rjk : k < fco} 
to all z G 5_E(77m, Tyivi), we obtain (6.48). □ 


It remains to prove Lemma 6.3. 
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Proof of Lemma 6.3. First we notice that E G I and (2.25) imply that for all sufficiently 
large N, the bounds (5.3)-(5.4) hold. Together with (6.41) they imply that 

^ 1 > ^ 1: (6-61) 


moreover, using (6.6) we also get 


1 


^ ^ 1- (6-62) 
Emniz) 

We start with (6.43). Thanks to symmetry, we only need to estimate \r\{z)\. By (6.61) 
we have 


\\G^{u;%{z))\\ = \\GA{io%{zm<i- 


Furthermore, UJgl'z) ~ 1 and Ima;g(S") > 1 imply TO^(a;^(z)) ^ 1. 
We next claim that 


E[tr(G^(w^(z))(A - z)Aa(z))] 


— 7V2^3 ’ 


(6.63) 


(6.64) 


for any e > 0 if > fVo(e) is large enough, uniformly for fj G [rirm ??m]- Assuming (6.64) and 
recalling the definition of i5^ and in (6.38)-(6.39), from (6.62) we get the first estimate 
in (6.43). 

Next, we prove (6.64). By the definitions in (6.37), we have 


E 


GiGj{uj%{z)){A - z)Aa{z)) 


= -E 
-E 


m[mH(z)] tT{G^{w%{z)){A - z)Gj^{z)) 
W[fBmHGj{Lo%(z))Gfiiz)) 


(6.65) 


We rewrite the two terms on the right side separately as covariances. 


E 


m[mH{z)] ti{G^{uj%{z)){A-z)Gf}{z)) 


respectively. 


E 


IE[/5(z)]tr {Gj{u;%{z))G^iz)) 


= Cov (mHiz),tT: (G^{uj%iz)){A - z)G^(z))) 
= GoY(^fg(z),tT{Gjiuj%{z))G^{z))^ , 


where Cov(X, Y) := E(]IE[W] • ]IE[F]), for arbitrary random variables X and Y. 

Given (6.42) and the uniform boundedness of mAasiz) from (5.16), we see that (6.12) 
is satisfied and we can apply Proposition 6.2 using different choices for Q. Together with 
Cauchy-Schwarz inequality |Cov(X, F)p < E|]IE[A']p •E|IIE[F]p, we get 


Cov (TOi/(z), tr (Gjiuj%{z)){A - z)Gg(z)^'j 
Cov (/b(?), tr (Gj(uj^(z))Gff(z))) 


1 

]Y2^3 ’ 
1 

N'^rp 


( 6 . 66 ) 


More specifically, for the first line of (6.66), we chose Q = I and Q = Ga(w^(z))(A — z); 
for the second line we chose Q = B and Q = Ga{^jJb^))^ where we also used the facts 
A = VAV* and B = UBU*. Here, we also implicitly used (6.63). Then, (6.64) follows 
from (6.66), which in turn proves (6.43). 

Next, using Proposition 4.1, (6.40) and (6.43) we immediately get (6.44). Moreover, 
since Itaoja^z) A"^’-^b{z) > Imz, we have \z — u]a{z) — ujBiz)\ > 1taujb{z) > 1. Together 
with (6.44) and (6.57) this implies (6.45). Notice that (6.42) together with the uniform bound 
on mAms imply the condition (6.12) in Proposition 6.2. Thus, finally (6.46) and (6.47) 
follow from (6.45) and the concentration inequality (6.13). This completes the proof of 
Lemma 6.3. □ 
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7. Two POINT MASS CASE 

In this section, we discuss stability properties of the free additive convolution /Xq, ffl fxp 
when both /Xq and fxp are convex combinations of two point masses. The analogous result 
to Theorem 2.5 is given in Proposition 7.2 below. Applications of that result in the spirit of 
Theorems 2.7 and 2.8 are then stated in Proposition 7.3 and Proposition 7.4. When we refer 
to the results in Sections 2-4, we will henceforth regard ni and ^2 as jia and respectively, 
unless specified otherwise. 

7.1. Stability in the two point mass case. Without loss of generality (up to shifting 
and scaling), we assume that 

Met = + (1 - , M/3 = C<5e + (1 - O'^o , ^*7^0, 

e,Ce(o,i], (7.1) 

Here we exclude the case = (~li \) since it is equivalent to = (1) \) 

under a shift. Note that the latter is a special case of Ma = M/S- 
Set 

£1 := min { ^ (l + ^ 

:= max { ^ (l + ^ 

4 := min { ^ (l + ^ + 

4 := max { ^ (^ + ^ + 
where we introduced 

r± :=e + C-2^C±\/4fC(l-0(l-C). (7-2) 

Note that i\ < ii < < 4- The following result, taken from [28], describes the regular 

bulk of in the setting of (7.1). Recall that denotes the density of (/x^ffl/x^)^'^. 

Lemma 7.1. Lei p,a and he as in (7.1). Then the regular bulk is given by 

=(^ 1 :^ 2 ) U (4,4), (7.3) 

in case /x„ ^ M/3, while in case jia = M/3 *■3 given by 

= (7.4) 

Proof. Choose the diagonal matrices A and B with spectral distribution /x^ = -I- (1 — 

^Ar)(5o and /xs = + (1 - Cn)4 respectively, with ■= XN\/N and (at := [CIVJ/A^, 

where [•] denotes the integer part. Recall from (7.1) that £, < C and + C ^ 1. From 
Theorem 1.1 of [28], we first observe that the 0 and 0 are eigenvalues of the matrix H = 
A+UBU*, U a Haar unitary, with multiplicities N{(n—£,n) and N{1—(n—^n), respectively. 
The remaining 2f^]\[N eigenvalues of H may be obtained via a two-fold transformation from 
the eigenvalues, {U), of a ^AfA^-dimensional Jacobi ensemble as 

rf:=i(l + 0±Y/(l-0)2+44), 


- 0)2 + 40r+) , i (1 + 0 - y'(l - 0)2 + 40r_) } , 

- 0)2 + 40r+) , i (1 + 0 - ^'(l - 0)2 + 40r_) } , 

- 0)2 + 40r+) , i (1 + 0 + a/( 1 - 0)2 + 40r_) } , 

- 0)2 + 40r+) , i (1 + 0 + \/(l - 0)2 + 40r_) } , 


j = 1,... XnN , 


(7.5) 
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and then identifying the eigenvalues of H as the set } U {tj } U {0, 9}. In addition, the 
weak limit of as iV —>■ oo, admits a density given by 


rf . _ 1 \/(r+ -a;)(x-r_) 

2 tt^ x{ 1 -x) 1 


X G 


(7.6) 


where and r_ are defined in (7.2). Since the limiting spectral distribution of H is 
given by ^q, ffl we see that (/Tq ffl agrees with the weak limit of the measure 

+ ^t“)> as — >■ oo. Using this information together with (7.5) and (7.6), one 

deduces that supp (^q, ffl = [^ 1 ,^ 2 ] U [^ 3 ,^ 4 ]. It then follows from the explicit form of 

the limiting distribution of the Jacobi ensemble that is bounded and strictly positive 

inside its support. This proves (7.3). 

In the special case fia = M/ 3 , we have ^2 = ^3 = 1 and thus supp (/Tq ffl fJ-p)^^ = [^ 1 ,^ 4 ], 
with ^1 = 1 — 2a/^(1 — () and £4 = 1 + 2-\/^(l — ^). In fact, the density of (/ia equals 

^ 1 \/(4-a;)(a:-£i) ^ ^ 


x[2 — x) ’ 
see (5.5) of [33] for instance. Then (7.4) follows directly 


□ 


Our main task in this section is to show the following result on the stability of the system 
(wq, W/ 3 , 2 ;) = 0 in the setting (7.1). Recall the definition of in (2.12). 

Proposition 7.2. Let fXa and /i /3 be as in (7.1). Let I C be a compact non-empty 

interval. Fix 0 < t/m < 00 . Then, there are constants k > 0, K < 00 and S < 00 , depending 
on the constants Y C, Vm and on the interval I, such that the subordination functions 
possess the following bounds: 


min lm.uja{z) > 2k, 
eSi(o,7)M) 

min Im. uja{z) > 2k, 
zeSxiO,r]M) 

(7.8) 

, ; M K 

max lwa(z)l<—, 
zeSi(0,r/M) 2 

, / M K 

.\ajpiz)\<—. 

z£Si{0,rfM) Z 

(7.9) 


Moreover, we have the following bounds: 

Y) Y Ta Y M/3; 

max V^j,^^^Y^a{z),w>p{z))<S. (7.10) 

^e5x(0,?7M) 

(a) Lf p,a = M/3, 

^Uoi,Ua{a}a{z)-:a:a{z)) < - — , ( 7 . 11 ) 

\z - 1| 

holds uniformly on iSz(0,?7m)- 

Remark 7.1. As an immediate consequence of Proposition 7.2 and (3.28), we obtain for 
Ma Y TP the bounds max^g^^/o,^^,) \uJa{z:)\ < 25, max,,g 5 ^(o,^^) \ujp{z)\ < 2S with I as 
in (7.10). For /i^ = pip, we get \uj'^{z)\ < uniformly on 5i(0,?7m) as in (7.11). 

Remark 7.2. In the case Ma = M/3, we note from Lemma 7.1 {c.f., (7.7)) that the point 
U = 1 is in the regular bulk However, m^^H^^(l + iO) is unstable under small 

perturbations. For instance, let 

TA = Ta = + (1 - C)<5o , MB = (^ - e)Ji + (1 - ,? + e)Jo , 

for some small e > 0. Then, according to Theorem 7.4 of [11], maHmb has a point mass eJi. 
Hence, even though (2.25) {i.e., di^{ptB, Tp) —0, as e —0) is satisfied, Bi^^H^^(z) contains 
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a singular part which blows up as 1^; — 1| = o{e). This explains, on a heuristic level, 

the bound in (7.11) and shows why the fj,a = fJ-p case at energy £1 = 1 is special even though 
the density is real analytic in a neighborhood of £1 = 1 . 

Remark 7.3. Consider a more general setup with /Xq, = ^da + 'fia and = (1 — ^)5b + pt/j, for 
some constants ^ € (0,1), a, 6 G R and for some Borel measures 'jla and Jip with /Xa(M) = 1—^ 
and Analogously to the discussion in Remark 7.2, we note that (a+ 6 +iO) 

is unstable under small perturbations. However, from Lemma 3.4, we know that the system 
{uia,u!p,z) = 0 is linearly 5-stable in the regular bulk under the assumptions of 
Theorem 2.5. That means, if neither p,a nor is supported at a single point and at least 
one of them is supported at more than two points, then the point E = a + h cannot lie in the 
regular bulk Thus, only in the special case p,a = M /3 with Ma as in (7.1), there is an 

unstable point, up to scaling and shifting given by £1 = 1 , inside the regular bulk 


Proof of Proposition 1.2. Estimates (7.8) and (7.9) follow from Lemma 3.2 and Lemma 3.3. 

To show statement (i), we recall from the proof of Lemma 3.4 that (wa, W/s, z) = 0 

is linearly 5-stable at {u!a,ujp) if 

|i-(^;jw/ 3 )-i)(^;,(‘^c)-i)i >c, (7.12) 

for some strictly positive constant c. We now show that (7.12) holds for the case Ma 7 ^ M /3 
in the setup of (7.1). Using henceforth the shorthand £b = £)i^, Fp = 5)^^, we compute 


Fc^iz) = 


zil-z) 


l-f-z’ 
Then it is easy to obtain 


Fp{z) = 


z{0 - z) 
e-ec-z’ 


FLiz) -1 = 


and 


\FLiz)-l\ = 


{1-f-zr’ 

lmFa{z) — Im 2 ; 
Im z 


F'p{z)-l = 


€ C+. 


g^(C-C^) 

{e-ec,-zf ’ 


(7.13) 


(7.14) 


Consequently, we have (c./., (3.18)) 


_ ImFgW-Imz 

^ Imz 


(Im 0 / 0 ( 2 :) — lm.z){\m.wp{z) — Imz) 


(7.15) 


Im o>a(z)Im UJp[z) 

for any z G C"*". Hence, for z G 5x(M0ibM) with some small but fixed ryo > 0 to be chosen 
below, we trivially get (7.12) from (7.15). It remains to discuss the regime z G 5 x( 0 ,mo)- 
Then (7.13) together with (2.5) implies that 


0 / 3(1 - 0 / 3 ) UJa( 0 -UJa) 


0/3(1 - O/3) 


= UJa + UJp — Z . 


1 — ^ — OJp 0 — 0( — (jJa ’ 1 — ^ — OJp 

Denote s := 1 — f — ujp and t '.= 0 — 0C, — oja- From (7.14) we then have 


(F'(o;3)-1)(A^(u;o)-1) = 


Using (7.16), some algebra reveals that 
1 1 f^+0-0C,-z 


st 






1 

st 


{sty 


1 0C + l-^-z 

-02{c-e)^ 0Hc-c)s ■ 


(7.16) 


(7.17) 


(7.18) 
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Owing to (7.15) and ~ (^C)^) > 0 (recall that C G (0,1/2] and 9 ^ 0), it 

suffices to show that 


|lm {st)\ > c, (7-19) 

in order to prove (7.12). Note that, from the definitions of s and t, together with (7.8) 
and (7.9), we have 

|lms|, |Imt| > c, Ishlil^C". (7.20) 

Since /ia there exists a positive constant d such that max{|^ — /■!, \9 — 1|} > d. It is 

then elementary to work out that 

max{|(^ - ^2) _ _ ^ 2^1 ^ 1 ( 2 ^ _ 20C + 0 - 1)1} > di, (7.21) 


for some positive constant c?i = c?i(C, Ci ^) > 0, since the special case {9, /) = (—1, 5 , ^) is 

also excluded in the setting (7.1). For brevity, we adopt the notation 

9C + l-^-z , ^ + 9-9(:-z 

9^c-e> ’ ic-e)t ■ 

Then, according to (7.18) we have 

11 11 , , 
Re - = Reil; - = Reip - g 2 (^_^ 2 ) ^ Im - = Im V'= Imv?. (7.22) 

If l(^ ~ ■C^) ~ ^^(C ~ C^)l > di holds in (7.21), then (7.22) implies that 

\Ke tj} — Re ip\ > d 2 , (7.23) 


for some positive constant d^ = d 2 (^,C,^)- For small enough rjo = ? 7 o(C, C,^*), then get 
{^ + 9 — 9(^ — E)Ret + O{r]o) (9(^ + 1 — ^ — E)Re s + O{r]o) 

- <PK - C^)kl»-' 

which, together with (7.20) and (7.23), implies that 

max{| 6 »C + l-C-F;|,|^ + 6 l- 6 lC-F;|} > ds , (7.24) 


for some positive constant da = dz{^, C, d). If, on the other hand, |(2^ — 20/ + 0 — 1)| > di 
holds in (7.21), we get (7.24) by triangle inequality. Either way, (7.24) follows from (7.21), 
for sufficiently small, but fixed, 770 > 0. 

Next, using (7.20) and (7.24), we see that there is a constant c > 0 such that, for 
sufficiently small 770 , for all z G 5x( 0 , 770 ), we have max{|Im(^|, |Im7/;|} > c. Since Ivcup = 
ImT/) by (7.22), (7.19) holds on 5i(0,77o). Therefore, (7.19) holds on all of 5 i(0,77m). So, if 
Ma M/3, the system (wc, W/j, z) = 0 is linearly S'-stable with some finite S. 

We next prove statement {ii) where /ia = /i,a and thus 0 = 1,/ = /. From (7.16), we see 
that oJa = w,a satisfies the equation 

^“(I ~ r,, . _ oc^ 

— 2u7a z. (7.25) 

-L s ^OL 

Solving (7.25) for a;a(z) we get 

Wa(z) = 07 / 3 ( 2 ;) = i (z - 1 + 2(1 - /) + v^(7^^T)2^^4//T^^) , (7.26) 


where the square root is chosen such that o;/ 3 (z) —>■ i\//(l — /), as z —/> 1. Substituting (7.26) 
into (7.17), together with the 0=1, / = /, s = t = l — / — o/a, yields 

4 (/-/ 2)2 


{F'^[ujp{z))-l){F'pMz))-l) = 


(z-l + V(z-l)2-4(/-a)^ 
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Then it is elementary to check that 

|l - (Ci(w/3(z)) - l){F'p{uJa{z)) - 1)1 > |z - 1|, z € S{0,rjM), 

which further implies (wa(z), u;/ 3 (z)) < \j\z— 1|. Hence (7.11) is proved. 


□ 


7.2. Applications of Proposition 7.2. Analogously to Theorem 2.5, we have two main 
applications of Proposition 7.2. The first one is the following modification of Theorem 2.7. 
Let \ia^ M /3 be as in (7.1) and let mAj IJ-b be arbitrary probability measures on K. Recall the 
domain Si{a,b) introduced in (2.13). For given (small) ir > 0, we set 

S^{a,b) := |z £ Si{a,b) : (;\z - 1| > max | A/dL(/iA, Ma). yTTb7T)}| ■ (7.27) 

Proposition 7.3. Let o,s in (7.1). Let X C be a compact non-empty 

interval. Let fiA, Lb be two probability measures on R. Fix 0 < tjm < oo. Then there are 
constants b > 0 and Z < oo such that the condition 


dL(MA,Ma) + dL(MBiM/3) < ^ (7.28) 

implies 

max |m^^H^3(z) - (z)| < Z (dL(MA,Ma) + dL(MB,M/3)) , (7.29) 

2e»Sz(o,77M) 

in case /Tq- ^ M/3? respectively 

^ \z-l\ (^L(MA,Ma) +dL(MB,Ma)) , (7.30) 

uniformly on 5^(0, mm) with o < Co, for some co > 0, in case pa = Lp- The constants b 
and Z depend only on the constants 9 and on the interval X, while co olso depends on b. 


Proof. Having established Proposition 7.2, the proof of (7.29) is the same as that of Theo¬ 
rem 2.7. To establish (7.30), we mimic the proof of Theorem 2.7 with S replaced by 
We only give a sketch here. Similarly to (5.9), using (7.8) and (7.11), we have with b in (7.28) 
sufficiently small that 

Tt,A,^,si.OJa{z),U)a{,z))< ^ ZG5J(0,7?m)- (7.31) 

|z - 1| 

As in the proof of Lemma 5.1, we rewrite the system ^^^^^^(u;aiz),u!a{z), z) = 0 as 
(oJa(z),uJa(z), z) = r(z) with ||r(z)|| satisfying the bound (5.8). From the uniqueness 
of the solution to (w^, ws, z) = 0 and (7.31), we get 

\w>a{z)-Wo,{z)\ < ||r(z)||/|z- 1|, 

|wb(z) - w„(z)| < ||r(z)||/|z-1|, zG5^(0,?7m), (7.32) 

via the Newton-Kantorovich theorem. Note that the inequality ||r(z)|| < — Ip is needed 

to guarantee that the first order term dominates over the higher order terms in the Taylor 
expansion of (ta^, ws, z) around (wq,, W/j, z). This is the reason why we restrict 

our discussion on the set 5^(0, ?7 m). In addition, thanks to (7.32) we see that (5.3) and (5.4) 
still hold with 5 i(0,mm) replaced by 5^(0 ,mm)- Then the remaining parts of the proof 
of (7.30) are the same as the counterparts in the proof of Theorem 2.7. □ 
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The second application of Proposition 7.2 gives the following local law for the Green 
function in the random matrix setup from Subsection 2.3.2. Fix any 7 > 0. We introduce a 
sub-domain of 5j(a, b) by setting 

5J(a, b) := S^ia, b) n [z € C : \z - 1\ > • (7-33) 

Proposition 7.4. Let be as in (7.1). Assume that the empirical eigenvalue distri¬ 

butions Pa, Pb of the sequences of matrices A, B satisfy (2.25). Fix any 0 < r]M < 00 , any 
small 7 > 0 and set . Let I C be a compact non-empty interval. Then 

we have the following conclusions. 

(*) If Pa ^ pp, then 

inax \mH{z)-mAmB{z)\< ./g/n ■ 

Z&SxiVmWM) Nrj'^IZ 

(a) If Pa = Pj3, then, for any fixed (small) c > 0, 

\mH{z) - mAfSBiz)] -< > 

uniformly on 5|(?7in,??M)- 

Proof of Proposition l.J^. Note that, in the proof of Theorem 2.8, the only place where we 
use the assumption that at least one of pa and pg is supported at more than two points 
is Lemma 3.4; in particular in (3.25). Hence, it suffices to mimic the proof of Theorem 2.8 
with Lemma 3.4 replaced by Proposition 7.2. Then the proof of the case pa ^ pp is exactly 
the same as that of Theorem 2.8. It suffices to discuss the case pa = pp below. 
Analogously to Corollary 5.2, with the aid of (7.31) and (7.32), we show that 

r^.4,MB(wA(2),WB(2;)) < , 2 ;g5|(0,7?m)- (7.34) 

Then, we use a continuity argument, based on Lemma 4.2 and Proposition 4.1 with S 
replaced by therein, to deduce from (6.40) that |a;f(z) — uii{z)\ < ||r^( 2 ;)||/| 2 ; — 1|, 

i = A,B, on S)^{r]ui, Pm)- The remaining parts of the proof are the same as in Theorem 2.8. 
This completes the proof of part {ii) of Proposition 7.4. □ 
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